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Abstract 

We study metastability and nucleation for the Blume-Capel model: 
a ferromagnetic nearest neighbour two-dimensional lattice system with 
spin variables taking values in {—1,0,-1-1}. We consider large but fi- 
nite volume, small fixed magnetic field h and chemical potential A in 
the limit of zero temperature; we analyze the first excursion from the 
metastable —1 configuration to the stable +1 configuration. We com- 
pute the asymptotic behaviour of the transition time and describe the 
typical tube of trajectories during the transition. We show that, unex- 
pectedly, the mechanism of transition changes abruptly when the line 
h = 2A is crossed. 

Keywords : Blume-Capel model, stochastic dynamics, metastability, nucle- 
ation. 
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Section 1. Introduction. 

Metastability is a relevant phenomenon for thermodynamic systems close to a first order 
phase transition. 

Let us start from a given pure equilibrium phase in a suitable region of the phase di- 
agram and change the thermodynamic parameters to values corresponding to a different 
equilibrium phase; then, in particular experimental situations, the system, instead of un- 
dergoing a phase transition, can still remain in a "wrong" equilibrium, far from the "true" 
one but actually close to what the equilibrium would be at the other side of the transi- 
tion. This apparent equilibrium, often called "metastable state" , persists untill an external 
perturbation or some spontaneous fluctuation leads the system to the stable equilibrium. 

For a general revue on metastability with particular attention to rigorous results see 
[PL1],[PL2]. 

There are strong arguments leading to the conclusion that neither metastability can be 
included in the scheme of Gibbsian formalism, which is confined to the description of the 
genuine stable equilibrium states (see [LR]); nor it can be directly described using extrap- 
olation beyond the condensation point, because of the presence of an essential singularity 
of the free energy (see the fundamental result due to Isakov [I]). 

Metastability is a genuine dynamical phenomenon. Its description on one side has a 
basic importance from the point of view of fundamental physics; on the other side it poses 
interesting new mathematical problems ([CGOV], [OS1], [OS2]). 

Since a general approach to non-equilibrium statistical mechanics is still missing, a cru- 
cial role is played by particular models of microscopic dynamics. It is remarkable that, 
quite recently, rigorous results have been deduced in this field by analyzing, in particu- 
lar, the geometry of the condensation nuclei as well as the possible coalescence between 
droplets. Notice that these aspects were totally absent in previous theories like the so 
called classical theory of nucleation (see [PL1]). 

In the recent years many progresses have been made in the understanding of the phe- 
nomenon of metastability in the framework of Glauber dynamics. By Glauber dynamics 
we mean a stochastic time evolution of a lattice spin system (in continuous or discrete 
time) whose elementary process is a single spin change and which is reversible (namely it 
satisfies the detailed balance condition) with respect to the Gibbs measure corresponding 
to a given hamiltonian. There is a certain freedom in chosing a particular dynamics sat- 
isfying the above mentioned requirements. A typical choice, that actually we will make in 
the present paper, is called "Metropolis dynamics" (see Eq. 2.6 below). 

The case of standard Ising model (spin +1 or —1, ferromagnetic nearest neighbour in- 
teraction), often referred to as Stochastic Ising model or Kinetic Ising model, has been 
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analyzed, in two dimensions, in [MOS] in connection to relaxation to equilibrium for arbi- 
trary large (and even infinite) systems close to the first order phase transition. 

A quite complete treatment appeared in [NS1] and [NS2] where J. Neves and R. Schon- 
mann analyzed, in the framework of the "pathwise approach to metastability" introduced 
some years ago in [CGOV], the phenomenon of nucleation for large but finite volume and 
small magnetic field in the zero temperature limit. 

In [SI] R. Schonmann, using an argument based on reversibility, described in detail the 
typical escape paths. 

Other asymptotic regimes, very interesting from a physical point of view and mathe- 
matically much more complicated, are considered in [S2], [SS]. 

In the same asymptotic regime as in [NS1], different hamiltonians have been considered 
in [KOI], [K02] where it has been shown that the typical path followed during the growth 
of the stable phase in general are not of Wulff type. Here by Wulff (shape) we mean 
equilibrium shape of a droplet at zero temperature namely the shape minimizing the surface 
energy for fixed volume. This non- Wulff growth seems to be an interesting phenomenon 
in the description of crystal growth. 

Let us now try to explain the nature of the mathematical difficulties related to our 
problem. We notice that in the above mentioned asymptotic regime the behaviour is 
similar to the one described by Freidlin and Wentzell in their analysis of small random 
perturbations of dynamical systems: the system typically performs random oscillations 
around the local minima of the energy and sometimes it goes against the drift following 
some preferential ways. In particular it is interesting to characterize the typical tube of 
trajectories during the first excursion from the metastable to the stable equilibrium. This 
first excursion can be seen as an escape from a sort of generalized basin Q of attraction of 
the metastable equilibrium. It turns out that many local equilibria are contained in Q and 
this more general situation goes beyond the approach developed in [FW] by Freidlin and 
Wentzell who were able to give a full description of the typical tube of escape only for the 
case of a domain D completely attracted by a unique stable point. 

In our more general case, as we will see, new interesting phenomena take place involving 
a sort of "temporal entropy". These phenomena are taken into account in [OS1], [OS2], 
where a complete description of the typical tube of escape is given in general. 

For attractive short range systems the main feature of the transition appears to be the 
formation of a critical nucleus with suitable shape and size. This critical droplet results 
from a competition between the bulk energy favouring the growth and the surface energy 
favouring the contraction. Only for large sizes and for particular shapes will the droplet 
tend to grow. 
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The present paper is devoted to the study of metastability and nucleation in the frame- 
work of a dynamical Blume-Capel model. It is a two-dimensional spin system where the 
single spin variable can take three possible values: —1, 0, +1. It was originally introduced 
to study the He 3 — He 4 phase transition. 

One can think of it as a system of particles with spin. The value a x = of the 
spin at the lattice site x will correspond to absence of particles (a vacancy), whereas the 
values a x = +1, —1 will correspond to the presence, at x, of a particle with spin +1, — 1, 
respectively. 

The formal hamiltonian is given by: 

H(a) = J (**-*v) 2 -AZ>l!-*Z>* > (1-1) 

<x,y> x x 

where A and h are two real parameters, having the meaning of the chemical potential 
and the external magnetic field, respectively; J is a real positive constant (ferromagnetic 
interaction) and < x,y > denotes a generic pair of nearest neighbour sites in Z 2 . 

In the following we will consider the system enclosed in a two-dimensional torus A. 
Let —1, and +1 denote the configurations with all the spins in A equal to — 1,0, +1, 
respectively. The structure of the set of ground states corresponding to different values 
of A and h will be discussed in Section 2. Now we only note that it is immediate to see 
that for A = h = the configurations —1, and +1 are the only ground states. It 
has been shown, using Pirogov-Sinai theory, that this phase transition persists at positive 
temperature T = 1//3 in the thermodinamic limit (see [B], [C], [BS] and [DM]). 

We will use as dynamics the Metropolis algorithm, in which the typical time needed to 
overcome an energy barrier AH is of order exp(/3AH). It will be defined in detail in the 
next section. 

We are interested to the case in which A and h are very small but fixed, the volume 
is large and fixed and T is very small; namely we move in the vicinities of the triple 
point h = A = T = 0. In particular we will consider the region h > A > where the 
the most interesting phenomena take place. The stable equilibrium, namely the absolute 
minimum of the energy, in this case, is +1 and we suppose to start with the system in the 
configuration — 1. We want to describe the first excursion between —1 and +1. It turns 
out that in the above region a direct interface between pluses and minuses is unstable and 
its appearence and resistance are very unlikely. 

The main effect which surprisingly and unexpectedly shows up is that two different 
mechanisms of transition between —1 and +1 take place for different values of the param- 
eters A, h. More precisely for < 2 A < h the transition takes place via the formation 
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of a suitable critical droplet of zeroes keeping growing untill it covers the whole volume. 
Subsequently, from the intermediate zero phase we have the nucleation of a critical droplet 
of plus spins driving eventually the system to the stable +1 phase. 

Conversely, forO<A</i<2A, the plus phase is created directly from the minus phase 
via the formation of a suitable critical nucleus, a sort of "picture frame" (see Fig. 3. 10), 
containing in its bulk the plus spins with a thin layer of zeroes separating the interior 
pluses from the sea of minuses. We want to stress that the line h = 2A, where this abrupt 
variation of the mechanism of nucleation takes place, has no meaning from the "static" 
point of view of the Gibbs states. The reason is that we are analyzing a region of the 
configuration space very unlikely at the equilibrium; but, on the other hand, this region 
and the form of the "energy landscape" (see Fig. 1.1) on it plays an important role in the 
relaxation from metastability. 



Fig. 1.1 Energy landscape of the Blume-Capel model. 

A first result that we obtain in the present paper refer to the computation of the 
asymptotic behaviour, for small temperatures, of the transition time (the life-time of the 
metastable state) . Then we pass to the characterization of the typical trajectories during 
the transition; we specify the geometrical sequence of droplets as well as the order of 
magnitude of the necessary time fluctuations during the growth of the critical nucleus 
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both for h < 2A and for h > 2A. To do this we exploit some general results contained in 
[OS1]. 

The model-dependent part of the work consists in the solution of a well specified se- 
quence of variational problems. The main difficulty is the determination of the "minimal 
global saddle" between —1 and +1 and of the set Q, the generalized basin of attraction of 
— 1. From this we will single out an optimal nucleation mechanism. We will analyze the 
energy landscape so precisely to exclude all the other possible mechanisms of transition. 
In particular we will show that any form of coalescence will be depressed in probability 
with respect to the optimal nucleation mechanism. 

The paper is organized in the following way: Section 2 contains definitions and results. 
In particular we state Theorem 1 concerning the asymptotics of the escape time. In Section 
3 we describe the local minima of the energy. In Section 4 we discuss supercriticality or 
subcriticality of droplets namely we determine their tendency to grow or shrink. In Section 
5 we prove a basic result on the height of different global saddles. In Section 6 we define 
the set Q and find the minima of the energy in its boundary dQ. In Section 7 we describe 
the typical tube of trajectories during the first excursion; then, using as preliminary results 
the propositions contained in the previous section we conclude the proof of Theorem 1; 
finally we state and prove Theorem 2 which refers to the typical tube. Section 8 contains 
the conclusions. Appendix 1 contains an explicit proof of a useful result about recurrence 
properties of a general class of Markov chains. 
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Section 2. Definitions and results. 

We start by describing the model that we want to study. The configuration space is 

Q A = {-1,0,+1} A , (2.1) 

where A = is a two-dimensional torus (a square with periodic boundary conditions) of 
side L. 

A configuration a is a function: 

a : A^ {-1,0, +1} . 
The energy associated to the configuration a is given by: 

H(a) = J (**-*») 2 -*5>2-ft5>* > ( 2 - 2 ) 

<x,y>CA xeA xeA 

where < x,y > denotes a generic pair of nearest neighbours sites in A and we suppose 
< A < h < J. We also introduce the following restriction 

2J 

A < 



2a 2 + a - 1 ' 



where a := j] the meaning of this condition will be clear later on (see (3.17)). 
The Gibbs measure in the torus A is given by: 

exp(-0H(v)) 

= , i 2 -*) 

where (3 represents the inverse temperature and Z\ is the normalization factor called 
partition function. 

We describe now the structure of the ground states corresponding to the different values 
of our parameters A and h. 

Let —1, and +1 denote the configurations with all the spins in A equal to —1, 0, +1, 
respectively. We have: 

for A = h = the ground state is three times 

degenerate, the configurations minimizing the 
energy are —1,0 and + 1 ; 
for h > and h > — A the ground state is + 1 ; 
for h < and h < A the ground state is — 1 ; 
for A < and A < h < — A the ground state is ; 
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Fig. 2.1 Ground states for the Blume-Capel model. 



for h = 0, A > : +1, -1 coexist. For h = A < : -1,0 coexist. For h = -A > : +1, 
coexist. These results are summarized in Fig. 2.1 where the coexistence lines are shown. 

We want now to introduce a dynamics in our model. It will be a discrete time Glauber 
dynamics namely a Markov chain with state space given by where 

1) the allowed transitions are between nearest neighbour configurations namely pairs £ 
and rj of configurations differing only in one spin: £ = rj x,b , with 



with b e {-1,0, +1}. 

2) It is reversible w.r.t. the Gibbs measure ji\ for the Blume-Capel model; namely the 
transition probabilities P(a, a') of the chain satisfy: 



Our choice is the so called Metropolis algorithm where the transition probabilities, for 
pairs of different configurations a, rj, are defined as 




Vy G A , y ^ x 
for y = x 



(2.4) 



ji A (a)P(a,a') = ji A (a')P(a' , a) . 



(2.5) 




a, n nearest neighbours 
otherwise 



(2.6) 
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where: 

+ f a if a > . . „ , n n , 

a + := < _ - . VaGR . (2.7) 

\ if a < v y 

The space of trajectories of the process is 

H:=(Ov) N . 
An element in S is denoted by u ; it is a function 

ui : N -> A . 

We often write cj = uo, <7i, . . . , cr f , . . .. 

We will call path an allowed trajectory namely: ui = o~o, a±, . . . , a t , . . . is a path iff aj and 
(Tj+i Vj are connected in the sense that (Tj+i = crj' 6 for some x E A and 6 G {—1, 0, +1}. 
We use the notation u> : a — > 77 to denote a path cj joining cr to 77. 

A path cj = o"o, o"i, . . . , a n is called downhill (uphill) iff if(a J+ i) < H(o~j) (H(aj+i) > 
H(o-j)) Vj = 0, 1, . . . , n — 1. We will use the convention that a downhill path can (and 
will) end only in a local minimum. 

A set Q of configurations, Q G Oa, is said to be connected iff for every pair of configu- 
rations a, 77 G Q, 3 a path u; : cr — > 77 such that uj C Q. 

We say that a configuration a is downhill connected to 77 iff there exists a downhill path 
uj : (J — ► 77. 

We will denote by M the set of all the locally stable configurations namely the set of all 
the local minima of the energy. More precisely: a G M iff for every x G A, b G { — 1, 0, +1} 
the corresponding increment in energy, given by 

A x , b H(a) :=H(a x > b )-H(a) (2.8) 

is positive. 

It is easy to see that in our model with the choice of the paremeters J, h, A that we have 

made, the quantity A Xj f ) H(a) will be always non-zero and this somehow simplifies some 
arguments. 

Given Q C Oa we define the (outer) boundary dQ of Q as the set: 

8Q := {a £ Q : 3a' G Q : P(</, a) > 0} , 

namely 

dQ := {a g" Q : 3x G A, 6 G {-1, 0, +1} such that a' = G Q} . (2.9) 
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We denote by U = U (Q) the set of all the minima of the energy on the boundary dQ of 
Q: 

U(Q) := {(EdQ: mm H(a) = H(()} (2.10) 

crEdQ 

and we define H(U(Q)) := if (f) with £ G U(Q). 

We denote by F = F(Q) the set of all minima of the energy on Q: 

F(Q) := {( E Q : minif (a) = H(()} . (2.11) 

crEQ 

Given a stable state a G M i.e. a local minimum for the energy, we define the following 
basins for a: 

i) the wide basin of attraction of a : 

B{a) := {( : 3 downhill path u : ( -> a} ; (2.12) 

ii) the basin of attraction of a given by: 

B(a) := {C : every downhill path starting from ( ends in a} , (2-13) 

B(a) can be seen as the usual basin of attraction of a with respect to the (3 = oo dynamics. 

iii) B(a) = the strict basin or attraction of a given by : 

B{a) := {C e B(a) : H{Q < H(U{B{a)))} . (2.14) 
We introduce now the useful the notion of cycle. A connected set A which satisfies: 

maxima) = H < min H(Q = H(U(A)) 

aeA CedA 

is called cycle. Notice that every local minimum for the energy is a (trivial) cycle. 

The following simple properties of the cycles are true. Their proof is immediate (see, 
for instance [OS1]). 

1. Given a state a G Oa and a real number c the set of all cr's connected to a by paths 
with energy always below c either coincides with Q\ or it is a cycle A with 

H{U{A)) > c . 



2. Given two cycles A\, A 2 , either i) A\ fl A 2 = or ii) A\ C A 2 or, viceversa, A 2 C A\ 
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We give now some more definitions: a cycle A for which there exists n* G U(A) downhill 
connected to some point a in A c , is called transient; given a transient cycle A the points 
r/* downhill connected to A c are called minimal saddles. The set of all the minimal saddles 
of a transient cycle A is denoted by S(A). 

A transient cycle A such that 3 a g A with H(a) < H(F(A)), there exists 77* G S(A) 
and a path oj : rf — > <7 6eZou> 77* (namely Vcr G cj : -H"(cr) < H(rj*)), is called metastable. 

For each pair of states a, ?? G Oa we define their minimal saddle «S(x, y) as the set of 
states corresponding to the solution of the following minimax problem: let, for any path u 

H(lo) := max .£/"(£), H ar) : = min H{ui) , 

find 

S{a, 77) := {( : H{() = H a , v ; 3 u : a ^ V , w 3 C, ^(w) = . 

One immediately verifies that a strict basin of attraction of a local minimum is a transient 
cycle. This case corresponds to a "one well" structure. More general cases involve the 
presence of "internal saddles" and correspond to a "several wells" situation. 

Given any set of configurations A C Oa, we use ta to denote the first hitting time to A: 

t a := inf{* > :a t E A} . (2.15) 

We use P v (-) to denote the probability distribution over the process starting at t = from 
the configuration n. 

We are interested in dynamics at very low temperatures. Namely, we will discuss the 
asymptotic behaviour, in the limit (3 — > 00, of typical paths of the first escape from — 1 to 
+1 for fixed h, A and A. 

Let us now better clarify the asymptotic regime in which we will operate: the volume 
|A|, the magnetic field h and the chemical potential A are fixed and we consider asymptotic 
estimates for (3 very large. This regime was studied in the case of standard Ising model 
in 2D by J. Neves and R. Schonmann in [NS1] where the point of view of the pathwise 
approach to metastability, introduced in [CGOV], was assumed. 

One can think, for instance, to take A very small, h = a A, a fixed positive number, |A| 
of order, say, of and (3 of order physically this corresponds to a regime in which, at 
the equilibrium, the energy dominates w.r.t. the entropy. 

In the above described situation the qualitative behaviour of our stochastic time evo- 
lution can be described as follows: the system will spend the majority of the time in the 
local minima of the energy. Sometimes it escapes from them but there is always a natural 
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tendency to follow a downhill path and an occasional, random and rather unprobable, 
uphill move. 

An important role will be played by the saddles separating different "basins of attrac- 
tion" (or generalized basins of attraction, see below) w.r.t. the (3 = oo dynamics. 

We will see that the local minima will correspond to particular geometric shapes that 
will be called plurirectangles (see Fig. 3. 8); we will analyze, in particular, the special saddles 
between "contiguous" local minima (see Lemma 5.1). 

A global saddle point is any configuration 

a eS(-l,+l) . 

In Section 6 we will see that the set of these global minima V are substantially different 
according to the values of the parameters A and h. 

1) For h < 2A we distinguish two cases: 

a) if 8 := P - 2J -^~V < %±A ? p i s Q f the form Vi, a given in Fig. 5.1 (the two critical 
lengths /* and M* are defined in (3.12) and (3.15)); namely it contains a "droplet" with 
external rectangle given by a square of side /* + 2; the internal shape given by a rectangle 
with sides /*, /* — 1, at a distance one from the external rectangle and with a unit square 
protuberance attached to the longest "free" side; the internal shape is full of pluses, the 
spins lying outside to the exterior rectangle are minuses; finally between the interior shape 
and the external rectangle there are zeroes. 

b) If 5 > 2jr, V is of the form Vxj, depicted in Fig. 5.1. is similar to V\, a but now 
the external rectangle has sides /* + 1, I* + 2 and internally we have a square with sides 
I* — 1 with a unit square protuberance attached to the shortest "free" side. 

2) For h > 2A, V is of the form Vi given in Fig. 5.1; namely it is given by a rectangle of 
sides M* and M* — 1, with a unit square protuberance attached to one of its longest sides, 
full of zeroes in a "sea" of minuses. 

We set: 

T := H(V)-H(-1) . (2.16) 

Let us now summarize our main results. 

We shall prove that the first excursion from — 1 to +1 typically passes through a con- 
figuration from V and the time needed for this to happen is of the order exp(/?r); this is 
the content of Theorem 1 that we are now going to state. 

Theorem 2 will characterize the typical trajectories of the first excursion. The proof of 
Theorems 1,2 and even the statement of Theorem 2 will need many more definitions and 
propositions. For this reasons they will be postponed to Section 7. 
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Theorem 1 is based, in particular, on Propositions 4.1, 4.2 and 4.3 given in Section 4. 
These propositions refer to the tendency of a given minimum r\ of the energy to evolve 
towards +1 or to —1 namely they establish under which conditions a droplet is supercritical 
or not. 

It will be crucial to introduce a sort of generalized basin of attraction of —1. Indeed we 
will reduce the proof of Theorem 1 to finding a certain set Q of configurations satisfying 
suitable properties. In order to explicitly construct this set Q we will need the results 
contained in Propositions 4.1, 4.2 and 4.3. This construction will be achieved in Section 
6. 

Theorem 1. 

Let f_i be the last instant in which at = —1 before r+i: 

f_i := max{t < r+i : a t = -1} . (2.17) 

Let 

f v := min{t > f_i : a t = V} ; (2.18) 

for every e > 0: 
i) 

lim P-i(f P < r+i) = 1 ; (2.19) 

ii) 

lim P_i(exp[/?(r -e)} < t +1 < exp[f3(F + e)]) = 1 . (2.20) 
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Section 3. Local minima of the hamiltonian H(a). 

In this section we want to analyze the geometrical structure of the local minima of the 
energy. 

For any configuration a G fl\ we denote by c + (cr), c~(o~) and c°(a) the union of all closed 
unit squares centered at sites x G A with a(x) respectivly equal to +1, —1 and 0. c + (cr), 
c~(a) and c°(a) decompose into maximal connected components c^' ' - , j = 1, . . . , k +,0, ~ . 

The centers of c^' ' form a ★-cluster in the sense of sites percolation, namely they are 
maximally connected components in the sense of the next nearest neighbours. The c^'°'~ 
will be simply called clusters. 

To any such c^' ' - we assign its rectangular envelope defined as the minimal closed 
rectangle R(cj ' ) containing it; if none of the rectangles R(cp°) is winding around the 
torus, we call the corresponding configuration acceptable. 

Let a be an acceptable configuration, we denote by 7^"'° the boundary of cp° Vj G 
{1, k + ' }; the internal component 7^'° of the boundary is defined as follows: let s be 
a unit segment of the dual lattice Z 2 + (|, |) belonging to 7j + '°, we say that s G 7*'° if 
and only if all the paths joining nearest neighbour sites of A and starting from the site 
adjacent to s and not in cp°, necessarily reach a site in cp° before touching the cluster cj 
winding around the torus. The external component 7^'° of the boundary of cp° is defined 
as 7*' \ 7j + '°- Of course 7*' can be empty. 

In order to construct the local minima of the hamiltonian we first prove that direct H — 
interfaces cannot exist in such configurations; in Fig. 3.1 we analyze the interaction of a 
minus spin with its neighbouring sites. We examine all the possible cases and we show 
that it is always possible to construct a lower energy configuration by changing the minus 
spin adjacent to the interface. 

Let a be an acceptable configuration such that there exists only one cluster of spins 
c° and no plus spins; it can be proved that 

{7 = 7 is a rectangle whose 
(3.1) 
sides are longer than two 

Indeed, if a is a local minimum and there exists a minus spin inside the cluster c°, then, as 
a consequence of the fact that c° does not wind around the torus, one has that necessarily 
it must exist at least one minus spin with at least two nearest neighbour sites occupied by 
spins (see Fig. 3. 2). 
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Fig.3.2 

This minus spin can be changed into + or by obtaining, in this way, a lower energy 
configuration, as shown in Fig. 3. 3; this is an absurd. 



u Jo -4J-(h-XJ r (/, A) 

,, I • 2// . \ ■ +8J-2h 



-0 

•---0- 
-0 



-2J-{h-X) 
+ +4J-2h 



Fig.3.3 

We can conclude that no minus spins can be inside c°, that is 7 = {0}. In a similar way 
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it can be proved that 7 is a rectangle and its sides are longer than two. 

The proof of the implication <= is in Fig. 3. 4, where it is shown that all the possible 
nearest neighbour configurations of a are at higher energy; in Fig. 3. 4 the modified spin is 
represented by a unit empty square. 
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Now let a be an acceptable configuration such that the following conditions are satisfied: 
there exists just one cluster c° of spins touching c~ , 7 is a rectangle, no minus spin is 
inside clusters of spins; all plus spins are in the cluster c + and 7" 1 " = 7 (see Fig. 3. 5). 
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With arguments similar to the ones used before, it can be proved that 



{<y+ = ^/+ i s a rectangle whose 
(3.2) 
sides are longer than two 



In the proof it is crucial that the energy of a configuration can be lowered by properly 
changing a spin having at most two zero spins and no minus spins among its nearest 
neighbour sites; all the possible situations are shown in Fig. 3. 6. 
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Fig.3.6 

Hence we have proved that configurations like the one in Fig. 3. 7 are local minima 
of H(a); these configurations are called birectangles and are denoted by the symbol 
R(L 1 ,L 2 ;M 1 ,M 2 ) where 



'Mi > Li + 2,M 2 > L 2 + 2 ifLi,L 2 >2 

Mi, M 2 > 2 if Li = L 2 = ' 1 } 
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It is easy to understand that the most general local minimum of H (a) is not a birectan- 
gle, but, rather, a more complicate configuration that we call family of plurirectangle (see 
Fig. 3. 8). It is an acceptable configuration satisfying the following conditions: 

i) there are k° clusters c°, c° o of spin touching c~; 

ii) 7°, •••,7fc are non interacting rectangles whose sides are longer than two; 

iii) in every cluster c°, there are kf clusters Ci~, •••,c~J" + of +1 spins; 

iv) Vj G {l,...,/c } 7^, 7 + fc + are non interacting rectangles whose sides are longer 

3 1 j 

than two. 

We have a single plurirectangle when k° = 1. 

We have used, above, the geometric notion of interacting rectangles: given two rectan- 
gles Ri and i?2 with boundaries on the dual lattice Z 2 + (|, |), we say that they interact 
if and only if one of the two following conditions occurs: 

i) their boundaries intersect; 

ii) there exists a unit square centered at some lattice site such that two of its edges are 
opposite and lie respectively on the boundaries of Ri and i?2- 

We have to compute the energy of such local minima as a first step in the description 
of their tendency to shrink or grow of the stables clusters. 
We say that a local minimum a is subcritical if and only if 

lim P a (T-i <r +i ) = 1 ; (3.4) 
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one of the main problems that we have to solve is to understand when a local minimum is 
subcritical. 

The energy of a birectangle R(Li, L 2 ; Mi, M 2 ) is 



H{R{L 1 ,L 2 ;M 1 ,M 2 ))-H{-1) 



(2Mi + 2M 2 ) J + (2Li + 2L 2 ) J - M x M 2 (h - A) - L x L 2 (h + A) . 



(3.5) 



The above formula can be easily generalized to the case of a general plurirectangle cr, char- 
acterized by the parameters Mij, M 2 j 7 e L 2 j^ Vj G {1, k } and Vz G {1, kj~} : 
with obvious meaning of the notation. One has 



H(a) - H(-l) 



J2{( 2M hJ + 2M 2,j)J - M 1J M 2J (h - A)+ 



k+ 



(3.6) 
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Now we consider a squared birectangle Q(L,M) := R(L, L; M, M), whose energy 
e(M, L) := H(Q(L, M)) - H(-l) is given by 



e(M, L) = AM J + ALJ - M 2 {h - A) - L 2 (h + A) 



(3.7) 



The graph of this function e : R 2 — > R is a paraboloid with elliptical section and downhill 
concavity, the coordinates of the vertex are 



M = 



2J 
h-X 



L = 



2J 

h + X ; 



(3.8) 



the level curves of e(M, L) are represented in Fig. 3. 9. 



2J/(h+A) 




2J/(h-X) 



Fig.3.9 

Let us consider a droplet Q(M, L) such that M < and L < if these conditions 

are satisfied e(M, L) is an increasing function of M and L, so we expect that this droplet 
will shrink. On the other hand if M > -^zj, since e(M, L) is a decreasing function of M, we 



expect that the external cluster of the droplet will grow; this suggests that M* 



2J 



+ 1 



h-\ 

is the critical dimension for the external cluster of a local minimum. After the growth of the 
external cluster, we look at what will happen to the internal one; with similar arguments 
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one can convince himself that L* 



2J 
h+\ 



+ 1 appears to play the role of the critical 



dimension. Obviously these two processes of growth cannot be inverted, in fact a plus spin 
droplet can "live" only inside a zero spin droplet. 

But it can also happen that the plus spin phase is reached directly, without passing 
through the zero spin phase; this happens if the droplet Q(M, L) grows moving along the 
line M = L + 2. In this case one can see that the system reaches the stable phase through 
a sequence of frames (picture frames). We call squared frame a birectangle C(l,l) := 
R(l, Z; / + 2, / + 2) with / > 2. The most general frame is a rectangular one (see Fig. 3. 10) 



C(h,l 2 ) :=R(h,l 2 ;h + 2,l 2 + 2) , 



where h,l 2 > 2. 



(3.9) 
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Fig.3.10 

Now we consider the energy of a squared frame e(l) := H(C(l,l)) — H(—l), using 
equality (3.7) we have 



e(Z) = -2hl 2 + Z[8J - 4(h - A)] + [8 J - 4(h - A)] ; 
the graph of this function is a concave parabola, whose vertex coordinate is 

2J- (h- A) 



(3.10) 



I 



h 



(3.11) 



We expect that C(Z, Z) will grow if Z > Z*, where 

"2 J — (h — X) 



V 



h 



+ 1 , 



(3.12) 
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otherwise it will shrink; hence /* should be the critical dimension of a squared frame. 

In order to describe the behaviour of a general birectangle R — R(L\, L2; Mi, M2), we 
must study growth and contraction mechanisms of a droplet; like in the Ising model these 
are mainly: growth of a (unit square) protuberance and corner erosion. But in Blume- 
Capel model the relevant local minima are made of two components, the internal and the 
external ones, and both of them can grow or shrink independently. The mechanisms of 
growth and contraction are explained in Fig. 3. 11, they corrispond to: 

1) creation of a + protuberance adjacent from the exterior to the internal rectangle; 

2) creation of a protuberance adjacent from the exterior to the external rectangle; 

3) erosion (H ► 0) of all but one + spin in a row or column of the internal rectangle; 

4) erosion (0 — > — ) of all but one spin in a row or column of the external rectangle. 

Their typical times are 

tl = e P[2J-(h+\)} f2 = e f3[2J-(h-\)] 

£3 = e P(h+X)(L-l) u = e f3(h-X)(M-l) ' ( 3 - 13 ) 

where L : = min{Li, L2} and M :— min{Mi, M2}. 
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Fig.3.11 

By comparing times t±, . . . , £4, we observe that the growth of an internal protuberance 
is always faster than the growth of an external one, indeed 

2J - (h + A) < 2J - (h - A) h<t 2 . (3.14) 
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Then we introduce the following critical dimensions 



(3.15) 



r 2J 1 r2j-2Ai 1 

M := +1 M:= — — +1 

h — A J /i — A 

whose meaning is explained below 

L < L* (h + A)(L — 1) < 2 J — (h + A): internal contraction is faster than growth, 

that is the internal component of the local minimum is (relatively) subcritical; 

L < L (h + A)(L — 1) < 2 J — (h — A): internal contraction is faster than external 

growth; 

M < M (h — \)(M — 1) < 2 J — (h + A): external contraction is faster than internal 
growth; 

M < M* <^> (h - A) (M - 1) < 2J - (h- X): external contraction is faster than growth, 
that is the external component of the local minimum is (relatively) subcritical. 



As we will see in the next section, another interesting length will be Iq 
We choose the parameters J, h and A in such a way that 



2J 

h+\- 



m + 1- 

2J+2A 2J 2J-2A 
h+X ' /i-A' /i-A 



2J ~fe l ~ A ' > an< ^ X are not i n teger, so that ambiguos situations, here and in the following, 
are avoided. 

The behaviour of our birectangle R depends on its dimensions, some of the possible 
cases are described below: 

L < L* and M < M*: both internal and external component are subcritical, R is 
subcritical; 

L < L* and M > M*: the internal component is subcritical, but not the external one, 
R is supercritical and the system starting from R will reach 0; 

L > L* and M > M*: both internal and external component are supercritical; R is 
supercritical and the system starting from R will reach +1 by passing through C(M\ — 
2, M.2 — 2) (internal growth is faster than external one); 

L > L* and M < M*: internal component is supercritical while external one is subcrit- 
ical, the future of the system starting from R depends on the relation M>M. 

Many different situations can take place, the last one is surely the most difficult but also 
the most interesting that we have to examine. 
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Growth and contraction of a frame are based on the same elementary mechanisms 
described before, but they take place in more than one step. The possible contrac- 
tion of a squared frame C(l,l) starts with the contraction of its internal component: 
our system typically first reaches the configuration S(l,l), increasing its energy of the 
quantity H(S(l,l)) — H{C{l,l)) = (h + X)(l — 1), and then the configuration R(IJ) := 
R(l-1, /; 1+2, 1+2), lowering its energy of the quantity H(S(l, l))-H(R(l, /)) = 2 J- (h+X) 
(see Fig.3.12). 
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Fig.3.12 

At this level it is not easy to describe the future evolution of the system: the internal 
component could continue to shrink or the external component could start its contraction; 
but we remark that the first step in the contraction of C(l, I) always involves bypassing of 
an energetical barrier whose height is (h + X)(l — 1). 

On the other hand the possible expansion of C(l, I) starts with the growth of an external 
protuberance: the system typically reaches the configuration G(l, I) by overcoming the 
energetical barrier H(G(l, /)) — H(C(l, I)) = 2 J — (h — X) and then it goes down to R(l + 
1, /) := R(l, /; / + 3, / + 2) lowering its energy of the quantity H(G(l, /)) - H(R(l + 1, /)) = 
(h — X)(l + 1) (see Fig. 3. 13). We have supposed, without loss of generality, that the growth 
is horizontal. 

As a consequence of the fact that it is always t\ < £2, the second step in the expansion 
of the droplet will be the growth of an internal protuberance: the system reaches the 
configuration S(l+1, 1) by overcoming the energetical barrier H(S(l+l, l))—H(R(l+l, /)) = 
2 J — (h + A) and then goes down to the frame C(l + 1, /) lowering its energy of the quantity 
H(S{1 + 1, /)) - H(C(l + 1, /)) = (h + X){l-l). 

In order to describe the future probable evolution of the system, starting from C(Z, /), 
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and establish its tendency to shrink or grow we have to distinguish the following four cases: 



L<1<1 



I <L => H(S(l,l)) < H(G(l,l)) 

(H(G(IJ))<H(S(IJ)) 
\ H{S(l,l)) < H(S{l + l,l)) 



I* <l, 1 + 2 < M 

r <i, m < 1 + 2 



J H(G(IJ)) < H(S(1 + 1J)) 
{H{S{l + l,l)) < H(S(IJ)) 

H(G(l,l)) < H(S(l,l)) 
H(S(l + 1, /)) - H(R(l + 1,1)) < H(G(l, /)) - H(R(l +1,1)) 



these four cases are illustrated in Fig. 3. 14. 

Even the analysis of growth and contraction mechanisms leads to conclude that I* is the 
critical dimension of a square frame. 

We close this section remarking that parameter A may be choosen sufficiently small, 
that is 

2 T 

A < — ^ (3.16) 

2a 2 + a - 1 

where a = j, so that the following inequalities are satisfied: 

1) L* + 1< r 

~~ ~ (3.17) 

2) L* < L < I* < I* + 3 < M < M* 



l//ebi)raio/2008 [25] 



3:25 



/ < L 



L<i<r 



S(l,l + 1)/ 




S(l,l+1) 



R(U + 1) 



C(l.l) 






I* < I, 1 + 2 < M 



I* < L M < I + 2 



Fig.3.14 
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Section 4. Supercriticality and subcriticality of local minima. 

In this section we want to prove rigorous results about supercriticality and subcriticality 
of local minima. Namely we want to give criteria to estabilish the natural tendency of the 
geometrical structures representing the minima for H to shrink or grow. We will first 
analyze the "frames" , then the generic birectangles and finally the plurirectangles. 

First of all we state the following proposition: 

Proposition 4.1. 

Let us consider the configuration C(li,l2); we set / := min{/i,/2} and m := maxjli,^}- 
Let e > 0, we have: 




where 




and 



m*(l) : 



2h 2J-(h-X) h + X 



I +1 . 



h — X h h — X 



Moreover 




where 




Finally 
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where 



ffl 2/ _ / e /9{[2J-(/i-A)]-(h-A)(m+l) + [2J 
T ±' ( £ ) : = | e /3[2J-(/ l -A)]±/3 £ 



otherwise 



Proof. 



Let us consider the frame C := C(/i, Z2) with / := min{/i, I2} < L, its basin of attraction 
B := S(C(/i,/ 2 )) and the relative boundary dB. Let us denote by Si G <9-B the set of 
configurations obtained by changing into zero / — 1 plus spin adjacent to one of the shortest 
sides of the internal rectangle of C (see Fig. 4.1); we claim that 

vam H(a) = H(Si) . (4.1) 

aEoB 



Remark. 

In the following we will consider: 

1. configurations a containing a unique droplet 7 with a given particular shape, size 
and location; for example a rectangle of zeroes (with given location and horizontal and 
vertical sizes) or a birectangle with given location and external, internal horizontal and 
vertical sizes. 

2. The equivalence class of all the configurations a' with a unique droplet 7' obtained 
from 7 by symmetries like rotations, translations, inversions w.r.t lattice axes and even 
displacements along sides of unit square protuberances. 

In the following, to avoid lengthy specification and to accelerate the exposition we often 
interchange the two above objects and we even use the same symbols to denote them. The 
reader will easily deduce the meaning of our statements from the context. 

For example sometimes we will denote by Si also a particular droplet obtained from a 
particular configuration in C by substituting one particular smaller internal side with a 
particular unit square protuberance. 

Let us now continue the proof of Proposition 4.1. 

In order to prove (4.1) we observe that, starting from C and considering all the possible 
uphill path, one is able to examine all the configurations in dB. The energy cost of 
all the possible first steps of the above mentioned paths are given in Fig. 4. 2; here we 
mark by a unitary square the site whose spin is changed and we denote by a couple 
of positive integer numbers the generic first step of our uphill path. We denote 
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by C^j the configuration reached after the step We observe that C 2}2 £ <9I? and 

that Ai^y > Ai7 2 2 V(z, j)^{(5, 1), (3,2), (2,2)}. Hence, all the paths whose first step is 
different from (5,1) and (3,2) lead to a boundary configuration whose energy is greater 
than H(C 2}2 ). 

Starting from C3 2 or C^i an uphill path can continue by following one of the ways shown 
in Fig. 4. 2 and in Fig. 4. 3. It can be easily shown that the steps (8, j) can be neglected as 
well. 
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In conclusion, only the paths made by steps (3,2) and (5,1) can lead to a configuration 
whose energy is lower than H(C 2 ,2)- 

Now, let a be an acceptable configuration such that the following conditions are satisfied: 
there exists just one cluster c° of spins which touches the sea of minuses namely the cluster 
c~ winding around the torus, no minus spins are inside c°; all plus spins are in a unique 
cluster c + included in c° and 7+ = 7 . If a E B then the following propositions are true: 

i) R(c°) = the external rectangle (l\ + 2) x (l 2 + 2) of the frame C; 
ii) R(c + ) = the internal rectangle h x l 2 of the frame C; 

Hi) the intersection of each one of the four sides of R(c°) with 7 contains at least a 
segment of length greater or equal to 2; 
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iv) the intersection of each one of the four sides of R(c + ) with 7+ contains at least a 
segment of length greater or equal to 2. 

We prove (i) by absurd: let us suppose that R(c°) is different from (Zi + 2) x (Z 2 + 2) and 
that 7 + is a rectangle. We can construct a downhill path which leads to a local minimum 
different from C by filling with spins the region R(c°) \ c+. Thus a & B, and this is an 
absurd, (ii) can be proved in a similar way. (Hi) is proved by absurd as well: suppose that 
the intersection between 7 and one of the sides of R(c ) contains only isolated intervals of 
length 1, namely there is a certain number of spins with three minus spins among their 
nearest neighbour sites. By changing this spins into —1 we construct a configuration at a 
lower energy level and characterized by a cluster of spins c'° such that R(c'°) is different 
from the rectangle (Zi + 2) x (Z 2 + 2); then there exists a downhill path which connects a 
to a local minimum different from C. Hence the absurd a B is obtained, (iv) is proved 
in a similar way. 

But, as we noticed before, all the uphill paths starting from C and leading to config- 
urations in dB with energy smaller than H(C 2,2) necessarily can only be made by steps 
(5,1) and (3,2). 

It is clear that, by virtue of the necessary conditions stated above, we cannot reach 
dB starting from C with less than I — 1 steps (5.1). On the other hand, since Si G dB, 
with more that I — 1 steps (5.1) we certainly get an energy larger that H(Si) and so a 
configuration which cannot be of minimal energy in dB. 

In this way we can only reach configurations with a unique cluster of pluses, so any 
boundary configuration with minimal energy is characterized by an external cluster c°, 
such that the intersection between 7 and all the sides of R(c°) is at least of length 2, and 
an internal cluster c + , such that the intersection between 7+ and one of the sides of R(c + ) 
has length 1 (see Fig. 4. 4). Among all these configurations it is easily seen that the one 
with lowest energy is 5i. 

In conclusion we have to compare H(Si) with if(C 2j2 ). Equality (4.1) follows from 
Z < L, H(Si) - H(C) = (h + X)(l- 1) and H(C 2 , 2 ) - H(C) = 2 J - (h - A). 

Now we want to apply to the description of the first escape from B the approach 
developed in [OS1], which is based on the properties of the above defined sets called 
cycles. 

It is easy to see that the basin of attraction B := B(C(li, Z 2 )) defined in (2.13) satisfies 
the following properties: 

i) B is connected; 
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Fig. 4. 4 

ii) Si C dB, and 

mill v H{a) = H(S 1 ), mm H(a) > H(S 1 ) 

aEoB /redB\Si 

Hi) Wr] G 5"i there exists a path u : i] — > C such that Vcr G u> \ {q} one has a G B and 
H(a) < H(S!). 

As it was noticed in [OS1] (see Proposition 3.4 therein), properties z), ii) and Hi) imply 
that the set B defined as the maximal connected set containing C and with energy less 
then H(Si) is a cycle with £1 belonging to its boundary dB. Moreover we notice, here, 
the following obvious properties: 

• any point r\ G dB necessarily is such that H{rf) > H(Si); 

• if ^(77) = H(Si) and r\ G" dB necessarily any downhill path starting from r\ ends in 
C. 

We recall that, given any set A C Oa, we have denoted by S(A) the possibly empty 
subset of U(A) (see (2.10)), which is downhill connected to A c ; S(A) was called the set of 
minimal saddles of A. We can write 

S(B) = S 1 . (4.2) 

From Proposition 3.7 in [OS1], from reversibility of the dynamics (see Lemma 1 in 
[KOI]) and from (4.2) we easily get that Va G B 

lim P a (a T 1 G 5i) = 1 . (4.3) 

p— >oo 

Since ^(^i) - H(C) = (h + X)(l- 1) we deduce that for every e > 

lim p c ^{h+m-D-Pe < Tqb < e p(h+X)(l-l)+0e ) = 1 (44) 
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Up to now we have described how the system reaches dB starting from C; now we want 
to describe its further evolution. 

Two things can happen: the system gets back to B or it goes to the birectangle R± := 
R(h — 1, 1; h + 2, 1 + 2) (see Fig. 4. 5); we have supposed, without loss of generality, that 



— 9AAAAAAAAAAAAAAAAAAAA — 9AAAAAAAAAAAAAAAAAAAA — OAAAAAAAAAAAAAAAAA^ 
-~r-0AAAAAAAAAAAAAA~r — r-OAAAAAAAAAAAAAA-T — r-6AAAAAAAAAAAAAA--T 

— i- 1 r i — r r r---r — r-r r— r 

— i- 1 r i — r r iiaa r — r-r iiaa r 

— i- 1 r— r-jA-r-r bae r a,-, i r r-r 

— I- 1 r i — i r 1 i — r-r hae i 

— i- 1 r i — i r 1 i — r-r r-r 

— r-nAAAAAAAAAAAAAS— r — r-IIAAAAAAAAAAAAAS— r — r-IIAAAAAAAAAAAAAS-T 

— nAAAAAAAAAAAAAAAAAAAS — IIAAAAAAAAAAAAAAAAAAAS — nAAAAAAAAAAAAAAAAAA 



In Appendix 1 we give a general argument showing that, with high probability, our 
process, possibly after many attempts, soon or later, will eventually get out of B U dB 
through Si reaching R\ before touching any other local minimum and: 



Now we have to describe the further evolution of our Markov chain starting from the 
birectangle R\. We denote by Bi := B(Ri) the basin of attraction of R\. Let us first 
consider the case min{/i — 1, /} = / (this is equivalent to suppose that C is not a squared 
frame). We denote by £2 the configuration obtained by changing into minus / + 1 of the 
spins on the "free" side of the external rectangle and by S3 the configuration obtained 
by changing into zero / — 1 of the plus spins of one of the shortest sides of the internal 
rectangle of R\ (see Fig. 4. 6). The following is true: 



Fig.4.5 



lim Pc{t Ri < 7+1) = 1 



lim P c (r Rl < e^+m-D+P^ = 1 



(4.5) 



min H(a) 

<TEdB 1 




(4.6) 
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Fig.4.6 

Equality (4.6) can be proved with arguments similar to those used in the case of the local 
minimum C(/i,/ 2 ) and observing that H(S 2 )-H(R 1 ) = (h-X)(l+l) and H(S 3 )-H(R 1 ) = 
(h + X)(l — 1) even though, in this case, there are other possible first steps (with high 
increment in energy). They are shown in Fig. 4. 7. 
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Fig.4.7 



With arguments similar to those used before we get that the typical time of first escape 
from dBi is of the order of e /3 [ mm<Teas i H ( (T )- H ( R i)] anc [ that the system hits for the first 
time the boundary dB 1 in S 3 if I < [j] + 1 and in S 2 if / > [j] + 1. Notice that if 
1 < j < 2 the integer [^] + 1 equals 2 so that £2 is preferred. 

We have that, with probability tending to 1 as f3 — > 00, our droplet continues its 
contraction: the system reaches another local minimum R 2 strictly contained in Ri, that 
is 

R 2 -< R 1 -< C ■ (4.7) 
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where we have introduced the following partial order relation in 

a -<n <^> a(x) < rj(x) Vx G A . (4.8) 

We also have that, given e > 0, e^ mmtre8B i fl"(o-)-H(.Ri)]+/3e j g an U pp er bound, in the limit 
(3 — > oo to the first hitting time to i?2 of the Markov chain starting from 

In conclusion we can say that the Markov chain starting from C visits smaller and 
smaller local minima untill it reaches the configuration — 1; this completes the proof of the 

statement Pc(t-i < r + i) /3 -^? ) 1. 

Each step of the shrinking process is characterized by a typical time tp whose asymptotic 
behaviour, exponentially in f3, is known in the sense that we control 

,. 1 i 

hm -\ogtp ■ 

p^oo p 

we say that: 

t\,t 2 Q are logarithmically equivalent <^ lim — = lim — ti . 

By using Markov property we can say that the typical time of the whole shrinking 
event is given by the largest time among all the partial shrinking times. Then the proof of 
Proposition 4.1 is completed in the case / < L when C is a rectangular frame. 

Next, we consider the case when C is a squared frame: the boundary configuration S3 is 
now the one represented in Fig. 4. 8, H(Ss) — H(Ri) = (h + X)(l — 2) and mm ae gB 1 H(a) = 
H(Ss) if / < [|f + |] + 1- We obtain results similar to those obtained in the previous 
case of a rectangular frame. 

9AAAAAAAAAAAAAAAA-- 
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1 1 IIAAAAAAAA I 

1 1 1 I 

1 1 1 I / + 2 

1 1 1 I 

1 1 1 1 

r-nAAAAAAAAAS'T- 

1 1 AAAAAAAAAAAAAAAE 

/ + 2 

Fig. 4. 8 
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Now we consider the frame C := C(l\, I2); we suppose that L < I := min{/i, I2} < I* and 
m := max{/i, I2} < m*(/); we denote by B the basin of attraction of the frame C and by 
dB its boundary. We denote by S 4 the set of configurations obtained by attaching a unit 
square protuberance (with a zero spin) to one of the four sides of the external rectangle of 
C (see Fig. 4. 9). By considering all the uphill paths starting from C it can be proved that 

min tf(a) = H(S 4 ), min H(a) > H(S 4 ) , (4.9) 

cr€oB aedB\S4 

namely 

U(B) = S 4 ; (4.10) 
we remark that H(S 4 ) - H(C) = 2 J - (h - A). 
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Fig. 4. 9 



Without loss of generality we suppose that / = I2 and m = l\. By arguments similar 
to those used before it can be proved that in a typical time e^ 2J ~^ h ~ x ^ the Markov chain 
starting from C, with high probability, will visit R2,± '■= R{h, h] h + 3, / 2 + 2) or i? 2 ,|| : = 
Rihih'ih + 2, 12 + 3). The symbol _L denotes the fact that the frame is growing in a 
direction perpendicular to its shortest side (see Fig. 4. 10). 
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Fig.4.10 

We denote by B 2t ± and B 2 ,\\ the basins of attraction of i?2,_L and R2,\\'i their boundaries 
are respectively denoted by 9i?2,j_ and dB 2 ,\\- By considering all the uphill paths starting 
from i?2,± and i? 2 ,|| we g e t that 



min H(a) = H(S 4 ) 

aedB 2 ,± 

min H(a) = H(S 4 ) 

creaB 2 ,|i 



(4.11) 



more precisely 

U(B 2 , 1 _) = S 4 , U{B U ) = S A . (4.12) 
Indeed the most relevant inequalities in the proof of (4.11) are the following ones 



(h + 1) >2J - (h - A) > 2J - (h + A) > {h - A)(Z + 1) 

(h + X)(l- 1) >2J - (h - A) > 2J - (h + A) > (h - A)(m + 1) 



(4.13) 



We remark that H(S 4 ) - H{R 2 ,i_) = (h-X)(l + l) and H(S 4 ) - H(R 2j \\) = (h-X)(m+l). 
In order to prove the first one of the equalities (4.13) we notice that 

I > L =>• (h + X)(l - 1) > 2J - (h - A) 
1 + 2 <M ^2J -(h + X) > (h- X)(l + 1) 
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In order to prove the second one we notice that 

/ > L =>■ m*(l) +2 < M 

and that 

m < m*(l) =>• m + 2 < M =>• (/i - A)(m + 1) < 2J - (/* + A) . 

Starting from -R2,_i_ or from i? 2 ,|| the system will typically go back to C before vis- 
iting other frames; these phenomena take place, respectively, in the two typical times 
e P{h-\){i+i) anc j e /3(/i-A)(m+i)_ j t appears clear that the system, before eventually leaving 
C to reach another frame, will wander, performing random oscillations, in the union of the 
basins B, B 2 ,± and B 2j \\. Then, in order to understand whether the frame will shrink or 
grow we have to describe its behaviour in a larger basin, containing B U B 2 ,± U B 2j \\ ■ This 
basin is denoted by V and it is defined as follows 

V := {q : every downhill path starting from r\ ends in C or R 2j ± or R 2j \\} ■ (4-14) 

We denote by S$ t ± and S^y the configurations obtained by attaching a unit square 
protuberance to the free side of the internal rectangle of R 2j ± and R 2i \\ (see Fig. 4. 11). By 
considering all the uphill paths starting from C, R 2j ± and R 2 ,\\, we are able to examine all 
the configurations in dV. We get: 

min H(a)=H(S 1 ), U(V) = S 1 

creoV 

The most relevant inequalities in the proof of equation (4.15) are 

/ < r (h + X)(l - 1) < [2 J -(h- A)] - (h - X)(l + 1) + [2 J -(h + A)] 
m < m*(l) ^(h + A)(Z - 1) < [2J - (h - A)] - (h - A)(m + 1) + [2 J - (/i + A)] ' 

(4.16) 

they mean, respectively, H(S5 j ±) > H(Si) and H(Ss ! \\) > H(Si). Of course it is always 
^(5 5 ,n) <H(S 5 ,±_). 

We notice that V is a sort of generalized basin of attraction of C; indeed it is easy to 
see that as a consequence of m < M* the "bottom" of V reduces to C in the sense that 
C are the only absolute minima of the energy in V and, as it is easy to see, starting from 
any initial configuration a G V our process, with high probability for large ft, will visit 
C before exiting from V. From V one can easily obtain, by suitably cutting in energy, a 
cycle having the same minimal saddles in its boundary: 
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(4.15) 
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Fig.4.11 



take the maximal connected set T> of configurations containing C with energy less than 
H(Si). Since it is easy to see that properties i), ii), Hi) of B = B(C(l\,l2)) are still 
verified with V in place of B for L < I < /*, m < m*(l) one immediately gets: S(T>) 3 Si; 
moreover, Vcr G V, 

lim P CT (cr T _i G Si) = 1 



and for every e > 



lim p c ( e 0(*+A)(I-l)-/3e < rw < e /3(^+A)G-l)+^ ) = 1 



We want to stress that the cycle T> is not the strict basin of attraction of any stable 
equilibium point but, rather, it has a several-well structure: it contains in its interior, 
beyond C, the equilibria i?2,± and B.2,\\', moreover it contains the internal saddles S4. The 
difference w.r.t the previous case of I < L (where we had to consider the cycle B(C) in 
place of V) is that now, not all the downhill paths emerging from a E V end in C and 
the system, before escaping from V will typically make many transitions, back and forth, 
between B(C) and B(R 2 ,±) B(R 2 ,\\) through S 4 . 

We consider, now, the frame C and suppose that L < I < I* and m*(l) < m < M — 2. 
We have H{Si) > H(S 5 ^) : H(S 4 ) < H(S 5 ^). With the usual arguments one can prove 
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that 

rmn H(a) = H(S 5 ^), U(V) = S^ ; (4.17) 
hence the frame C is supercritical and the system starting from C will hit +1 in a typical 

time e /H[2J-(h-A)]-(ft-A)(m+l)+[2J-(h+A)]}_ 

In the case L < I < I* and m*(Z) < M — 2 < m it can be proved that mino-ggp = 
H(S 5j \\), U(V) = S 5j \\ and H(S 5j \\) < i/(S , 4 ). Hence the frame is supercritical, but the 
typical escape time is e 0l 2 J-(h-X)] . We remark that in this case H{Sa) < H{S{) < H(S 5j ±) 
and H(S 5j \\) < H(S 4 ) (see Fig.4.12). Hence the set 

V:={<j e V- H(a) < H(S 4 )} 

is a generalized cycle like the set Qi defined in [OS1]. V is a set of cycles communicating 
through the minimal saddles in S4. Starting from V the system, before leaving V will not 
necessarily visit all the cycles contained in V with energy less than H{S 4 ). 

The system can leave B either through S^± or through S 4 ^. In the first case it will 
enter into B 2 ± visiting all the configurations of the cycle: 

B 2t± := {a G B 2 ,± : H(a) < H(S 4 )} 

before leaving B 2j ± and passing again through S 4j ±. In the second case it will directly get 
out of V. 

In the case /* < / < M — 2 and to < M — 2 we have that min cre gx? H (a) = H(V 5 ^\), 
U(V) = S 5j \\ and H(S 5 ^) > H(S4). Hence the frame is supercritical and the typical escape 
time is e / 3 {[2«/-0-A)]-(/i-A)(m+i)+[2j-(h+A)]} _ j n cage mos t important inequalities 

are 

(h + X)(l- 1) >[2J -(h- A)] - (h - X)(l + 1) + [2 J - (* + A)] > 

(4.18) 

>[2J -{h- A)] - (/i - A) (to + 1) + [2 J -{h + A)] 

we remark that #(Si) - #(C) = (/i + A)(Z - 1), #(S 5>± ) - H{C) = [2 J - (h - A)] - (h - 
X)(l + l) + [2J-(h + X)} and H(S 5j]l )-H(C) = [2J -(h-X)]-(h-X)(m + l) + [2J -(h + X)}. 

In the case /* < I < M — 2 and M — 2 < to we have that mm ae Qx> H(a) = H(S 5 ^) 
and H(S 5 ^) < #(64). Hence the frame is supercritical and the typical escape time is 
e /3[2j-(h-A)] j n cage we k ave that j) con tains again a generalized cycle. 

We remark that in the supercritical cases discussed above, namely for / < M — 2, the 
growth of a rectangular frame is asymmetric. The frame grows in a direction parallel to 
its shortest side towards a squared frame. Notice that the same tendency to be attracted 
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Fig.4.12 

by a squared shape is present also in the contraction of a subcritical frame which, as we 
have seen above, prefers to shrink in the direction orthogonal to its smallest side. 

Finally we consider the case I > M — 2. It can be proved that H(S^^) < .ff (S^jJ < 
H^Sa) < H(Si), hence the frame is supercritical and the typical time is e^ 2J ~^ h ~ x ^ . In 
this case the growth process is symmetric, similary to what happens in the stochastic Ising 
model for any supercritical rectangle. This concludes the proof of Proposition 4.1. □ 

Remark. 

In the following to avoid lengthy and repetitions we will often use short expressions 
like: the external rectangle shrinks in a direction perpendicular to its shortest sides instead 
of: by a comparative analysis of the possible barriers of energy, namely looking at the set 
of minimal saddles of a suitable (possibly generalized) basin of attraction, we know that 
with a probability tending to one as (3 tends to infinity the external rectangle shrinks in a 
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direction perpendicular to its shortest sides. 



In Proposition 4.1 we have stated conditions of subcriticality and supercriticality for 
frames, now we state similar results for birectangles 

Proposition 4.2. 

Let us consider a birectangle R := R(Li, L2; Mi, M2), let L := min{Li,L2}, L : = 
max{Li,L2}, M := min{Mi,M2} and M := max{Mi,M2}. If one of the conditions 
1-5 is satisfied 

1) L < L*, M < M* 

2) L > L*, M <M, L + 2 < M, L <L 

3) L> L*, M < M, L + 2 < M, L < L < /*, L < m*(L) 

4) L > L*, M <M, L + 2 > M, M - 2 < L 

5) L > L*, M <M, L + 2 > M, L<M-2<l\L< m*(M - 2) 

then 

lim Pr(t-i < r+i) = 1 . 
If one of the conditions 6-11 is satisfied 

6) L > L*, M <M < M* 

7) M > M* 

8) L > L*, M <M, L + 2 < M, L > I* 

9) L > L*, M < M, L + 2 < M, L < L < l\ L> m*(L) 

10) L > L*, M < M, L + 2 > M, M — 2 > /* 

11) L > L*, M < M, L + 2 > M, L < M — 2 < /*, L > m*(M — 2) 

then 

lim Pr(t+i < r_i) = 1 . 

I ':/,/,[. 

Without loss of generality we can assume M = Mi. Let us denote by B := B(R) the 
basin of attraction of R and by dB its boundary; first af all we have to find the minimum 

1/ febbraio/2008 [42] 4:42 



of the energy on the boundary dB. We examine all the uphill paths starting from R, but 
the relevant ones are those made of steps of the kinds (2, 2), (3, 2), (5, 1) and (10, 1) (see 
Fig. 4. 2 and Fig. 4. 7). The boundary configurations Si, S%, S3 and S4 reached by the uphill 
paths described above are represented in Fig. 4. 13, togheter with the energy differences 
AHi = H(Si) — H(R) Vz = 1, ... ,4. We remark that certainly if both the shortest sides 
of the external rectangle are not "free" , then at least one of the longest sides will be free; 
in this case AH 2 = (M — l)(h — A). 

Now, we consider the case (1): as a consequence of the subcriticality of the internal and 
the external rectangle one has 



But we cannot say anything about the inequality (L — l)(h + A) > (M — l)(h — A), without 
specifying more conditions on L and M; therefore the minimum of the energy on dB is 
given either by S2 or S3, depending on the values of L and M. Thus, the birectangle R 
is always subcritical and we can express the typical time needed by the system starting 
from R to hit —1 as max{e^ i_1 ^' 1+A ^ e 0(. M - 1 )( f >—X)y^ Similar results are obtained if one 
supposes that the shortest side of the external rectangle is not "free" . 

Now, we suppose L > L*, M < M and L + 2 < M (see Fig. 4. 14): the internal rectangle 
is supercritical, namely 2 J — (h + A) < (L — l)(h + A), and the external one is subcritical, 
namely (M — l)(h — A) < 2 J — (h — A); moreover 



Then the minimum of the energy on the boundary dB is £2 . The external rectangle shrinks 
in a direction perpendicular to its shortest sides untill it becomes a squared rectangle, then 
the shrinking process goes on in both directions untill the frame C(Li, L2) is reached in a 
typical time e P(M-i)(h-\) _ 

Even in the case L > L*, M < M and M < L + 2 < M (the longest sides of the internal 
and the external rectangle are necessary parallel) the system, starting from R, reaches the 
frame C(Li,L2). Indeed the external rectangle shrinks along the direction perpendicular 
to its shortest sides untill this process is stopped by the internal rectangle (see Fig. 4. 15). 
In other words this appears when the configuration R(L, L; M, L + 2) is reached (we have 
supposed, without loss of generality, that L\ = L). At this point the external rectangle will 
begin to shrink along the direction perpendicular to its longest sides, because L + 2 < M 



L<L* =^ (L- 
M < M* ^ (M 



l)(h + A) < 2J - (h + A) < 2J - (h - A) 
-l)(h-X) <2J-(h-X) 



(4.19) 



M < M =>■ (M - l)(h - A) < 2J - (h + A) . 



(4.20) 
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Fig.4.13 

and then (L + l)(h — A) < 2 J — (h + A). Hence, the system, starting from R, reaches the 
frame C(L U L 2 ) in a typical time max{e^ +1 ^- A \ e P(M-i)(h-X)y = e /3(L+i)(/ l -A)_ 

Then, we can conclude that in the cases (2) and (3) the frame C(Li,L 2 ) is eventually 
reached, but C(L 1 ,L 2 ) is subcritical, hence R is subcritcal, as well. For similar reasons in 
the cases (8) and (9) the birectangle R is supercritical. 
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Fig.4.15 

The typical shrinking time is given by 

max{e' 3 ( M - 1 )( fc - A )e« L - 1 )(' i+A )} if L + 2 < M 
m^{e^ L+1 ^ h - x \e^ L - 1 ^ h+ ^} if L + 2 > M ' 

With similar arguments it can be shown that in cases (4) and (5) the system, starting 
from R, hits C(M— 2, L) in a typical time e 2J ~( h + x )_ Hence, the birectangle R is subcritical 



l//ebi)raio/2008 [45] 



4:45 



and the typical shrinking time is e 2J -(M-A)_ i n the cases (10) and (11) the birectangle R 
is supercritical, as a consequence of the supercriticality of the frame C(M — 2, L). 

With arguments similar to those used before it can also be seen that in the case (6) the 
birectangle is supercritical since it first evolves towards the frame C(M — 2, M — 2) which 
is a supercritical frame since M — 2 > /* with our choice of the paremeters. 

Finally, in the case (7), the birectangle is easily seen to be supercritical. Indeed it 
follows from an argument similar to the corresponding one valid for the standard Ising 
model that starting from a configuration with M > M*, we get before +1 in a time of 
order e ^ 2J ~^ h ~ x ^ with high probability for large (3. Then, starting from we tipically 
follow an Ising-like nucleation path (see [NS1], [SI]) leading to +1 through the saddles 
5(0, +1). These saddles are given by configurations with precisely one cluster of pluses 
(in the sea of zeroes), this cluster is given by a rectangle L* x (L* — 1) with a unit square 
protuberance attached to one of its longest sides. It is immediate to verify that 

H(S(0,+1))<H(V) . 

The proof of Proposition 4.2 is complete. □ 

We consider, now, a plurirectangle R. We denote by Mi and M 2 the lengths of the sides 
of the external rectangle, by L\^ and L 2 ,i Vz = 1, k + the lengths of the sides of the k + 
internal rectangles Rf and we define M := min{Mi,M2} and Lj := minjLi^, £2,1} = 
1, k + . In order to state conditions of subcriticality and supercriticality for such configu- 
rations, we must introduce the rectangle R + defined as the rectangular envelpe of the union 
of all the internal supercritical rectangles. We denote by L 1R + and L 2R + the lenghts of 
its sides and we define L R + := mm{L 1R + , L 2j -R+} an d L R + :— max-fl^ jR +, L 2yR +}- Sup- 
pose that 3i E {1, 2, k + } such that Li > L*, we denote by R the birectangle obtained 
by removing all the internal rectangles and by filling up with plus spin the rectangle R + . 
Finally we state the following proposition 

Proposition 4.3. 

If one of the two following conditiones is satisfied 

1) Li < /*Vz = 1, k + and M < M*; 

2) 3i G {1, 2, k + } such that Lj > L* and R is subcritical; 
then 

lim Pr(t_i < r+i) = 1 . 

Proof. 
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Let us consider the case (1): we prove Proposition 4.3 by describing the shrinking 
process. 

First af all the internal rectangles whose sides are such that (Lj — l)(h + A) < (M — 
l)(h— A) shrink in a typical time e^^ Li ~ 1 ^ h+ ^ . We denote by the rectangular envelope 
of the union of all the "surviving" rectangles R^Vi G 1^ C {1, k + } and by l/ 1 ) its 
longest side. 

At this point the external rectangle starts shrinking (if it can). If ZA 1 ) < M — 2 this 
contraction ends when the external rectangle reaches R^ (see Fig. 4. 16). 
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Fig.4.16 



Let us define Z min := min ieJ (i) {Li}: the internal rectangle Rf such that Li = L min 
starts shrinking and loses a slide of lenght Li = L min . There are two possible situations 
(see Fig. 4. 17): after this contraction the external rectangle has a "free" side or not. In the 
first case the external rectangle loses another slice and a configuration of the type described 
in Fig.4.16 is reached. In the second case the internal rectangle goes on shrinking untill it 
disappears, and a configuration like the one in Fig.4.16 is reached, as well. In both cases 
the plurirectangle goes on shrinking by the mechanism described before untill it disappears, 
hence in the case l/ 1 ) < M — 2 the plurirectangle R is subcritical. 

Now, we consider the case L^ > M — 2. During the second phase of the contraction 
the system reaches a configuration characterized by an external rectangle whose sides 
are M and L^ + 2. The "free" side of the external rectangle is eventually L^ + 2. 
If (LW + l)(h - A) < (Li - l)(h + A) Vz G the external rectangle shrinks in a 
direction perpendicular to its "free" side untill it reaches R^ 1 '; and then the shrinkig 
goes on as we have described before. If there exists an internal rectangle R~l such that 
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Fig.4.17 

(L^ 1 ' + l)(h — A) > (Li — l)(h + A) it disappers before anything else can happen. Then 
the contraction goes on as described before. In conclusion we have proved that in the case 
(1) the plurirectangle R is subcritical. 

In the case (2) the proof of Proposition 4.3 can be achieved with arguments similar to 
those used in the case (1). □ 
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Section 5. Comparison between special saddles. 

Let us consider a subcritical frame or birectangle; we say that such a configuration is 
almost-supercritical iff it can be transformed into a supercritical minimum by attaching 
to one of its internal or external sides a whole slice. By attaching a slice to an internal or 
external side of a birectangle (or, in particular, of a frame) we mean transforming from — 1 
to the value of the spins in the row or column adjacent externally to this side. "Removing 
a slice" is the inverse operation of "attaching a slice" . 

Let us consider, now, a supercritical frame or birectangle; we say that such a configura- 
tion is just-supercritical iff it can be transformed into a subcritical minimum by removing 
a whole slice from one of its internal or external sides. 

Let us consider an almost supercritical frame or birectangle, we denote by u the inter- 
nal or external side such that by attaching to it a whole slice we obtain a supercritical 
configuration. We call special saddle the configuration obtained by attaching to u a plus 
unit protuberance, if u is an internal side, or a zero unit protuberance, if u is an external 
one. 
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Fig.5.1 

Let us consider the set V := (V\ UP 2 ) C Oa with V\ and P 2 the set of special saddles 
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shown in Fig. 5.1, where we have used the following definition 



i:=/--?^£z*>e]0.1[. (5.1) 
We state the following lemma: 

Lemma 5.1. 

For any special saddle S £ V it there exists S* E V such that 

H(S) > H(S*) . 

Before starting the proof, we observe that the frame C(l*, I*) is supercritical and C(l* — 
1,1* — 1) is subcritical for any choice of the parameters A and h; indeed it can be proved 
that 

m*(l* — 1) > I* for any value of h and A, (5.2) 

(see (5.5)). On the other hand we remark that the criticality of the frame C(l* — 1,1*) 
depends on the value of the real number S defined in (5.1). By comparing the energies of 
the two configurations shown in Fig. 5. 2 one can easily convince himself that 

C(l* -1,1*) subcritical iff 5 < + A 



C(l* — 1, I*) supercritical iff 5 > 



2h 
h + X 
2h 



we observe that e]0, 1[ if j > 1. This explains the reason of the twofold definition of 
the configuration V\. 
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Fig.5.2 
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Proof of lemma 5.1. 

Let us suppose that 5 < One can prove that for any / such that L </</* — 1 

m*(l) > 1* + 1 . (5.3) 

First of all we observe that m*(l) is a decreasing function of I, more precisily one can easily 
prove that 

m*(l - 1) > m*(Z) + 1 VZ G [Z, r - 1] . (5.4) 

Therefore in order to get a lower bound to m*(Z) it is sufficient to evaluate m*(l* — 1); 
with some algebra one can easily obtain 



m *(l* _!) = /* + 

Then, 



(i-*) 2! ' 



/i - A 



(5.5) 



ff <^ = ,(i_^>i^ ro -(r-i)>r + i ; 

this completes the proof of inequality (5.3). We remark that the validity of the equations 
(5.4) and (5.5) does not depend on the value of the real number 5. 

Now, in order to prove Lemma 5.1 we have to examine all the possible special saddles. 

Case CI. 

We consider the special saddle C\(m) in Fig. 5. 3. 
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Fig.5.3 
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It can be easily shown that H(Ci(m)) is an increasing function of m G [/*, m*(l* — 1) — 1], 
indeed H(Ci(m + l))-H(Ci(m)) = (h + X)-2hd > by virtue of the hypothesis 5 < 
Hence, 

H(d(m)) >Vi Vm G [l*,m*(l* - 1) - 1] ; (5.6) 
we observe that the equality is verified in (5.6) iff m = /*, that is Ci(m) = V\. 
Case C2. 

We consider the special saddles C2, a (l) an d C2,b(0 in Fig. 5. 4. We remark that the 
configuration obtained from C2,b(l) by removing the protuberance is subcritical because 
m*(l - 1) > m*(/) + 1. 
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Fig.5.4 



We have that H(C 2 , a (l)) < H(C 2 , b (l)); indeed H{C 2 , a {l)) - H(C 2 , b (l)) = {h+ A)(Z - 
m*{l)) and / — m*{l) < 0; the last inequality is a consequence of the fact that / < /* and 
of equation (5.2): m*(Z) > m*(l* — 1) > I* > I. 
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We observe that H{C 2 ^ a {l)) is a decreasing function of I: 

H(C 2 , a (l + 1)) < H(C 2 , a (l)) VZ E [L, /* - 2] . (5.7) 

Indeed, it is not difficult to show that H(C 2 , a (l + 1)) - H(C 2 , a (l)) = (h + A) + 
2h (I* -I- 5) (m*(l + l)-m*(l))-2h (m*(l + 1) - I* + 5) and by observing that I* -1-6 > 
+1, m*(Z+l)-m*(Z) < —I andm*(/ + f)-/*+5 < we obtain H{C 2 , a {l+l))-H{C 2 , a {l)) < 
(h + A) — 2h = A — h < 0. This completes the proof of the inequality (5.7). 

Since H(C 2 ^ a (l)) is a decreasing function, we have to compare the energy of the two 
configurations C 2ja (l* — 1) and Pi; by a direct calculation one obtains H(C 2 ^ a (l* — 1)) > 
H(V!). 

Case Bl. 

We consider the special saddles Bi jU (M, M;L) and #i 5 (,(M; L, L) in Fig. 5. 5 (here and 
in the following we use the notation introduced in Proposition 4.2 to label the internal 
and external sides; we use L, L, M and M to denote the dimensions of the birectangle 
obtained by removing the unit protuberance of the special saddle). 
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(a) (b) 

Fig. 5. 5 (a) The internal horizontal dimension is L* — 1 and the vertical one, L, is such that 
L>L*. The external vertical dimension is M and the horizontal one M is such that M<M<M* . 
If we choose the parameters h and A such that M=M*, then the special saddle (a) does not 
exist, (b) The external horizontal dimension is M— 1 and the vertical one is M. The internal 
dimensions L and L are such that L > L* and by removing the external unit protuberance one 
obtains a subcritical birectangle. 

We observe that 

H(B ha (M, M;L))> H(B ha (M, M; L*)) (5.8) 
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for every possible choice of the positive integer numbers M, M and L. This is an obvious 
consequence of the fact that L = L* — 1 < L* and M < M*. 

Now, we transform £>i )(1 (M, M; L*) into V\ in several steps and we evaluate the energy 
cost AHi of each step. 

• B lt a(M,M;L*) -> R(L* - 1,L*;M,M), A#i := H(R(L* - 1,L*;M,M)) - 
H{B lia {M, M-L*)) = -[2J -{h + A)]; 

• #(L* - 1, L*; M, M) -> #(L* - 1, L*; Z* + 2, Z* + 2), AH 2 < because the external 
rectangle is subcritical and M > I* + 2 (see inequalities (3.17)); 

• R(L* - 1, L*; Z* + 2, Z* + 2) -> R(L*, L*; Z* + 2, I* + 2), Aif 3 < because a whole 
internal slice of lenght L* has been attached to the internal (relatively) supercritical 
rectangle; 

• R(L*, L*; I* + 2, Z* + 2) -> i?(Z* - 1, /*; /* + 2, Z* + 2), A# 4 < because the internal 
rectangle is supercritical and L* < I* — 1; 

• #(/* - 1, /*; I* + 2, I* + 2) -> V u AH 5 = 2J - (h + A). 

One has that J]J , =1 A^ < 0, hence H(B lja (M, M; L*)) > if (Pi). This inequality and 
(5.8) lead us to the conclusion that 

H(B 1 , a (M,M;L))>V 1 (5.9) 

for every possible choice of M, M and L. 

In order to carachterize the special saddle Bi^{M\ L, L), we have to distinguish two 
possible cases. 

Case (i) L + 2 > M: the birectangle R(L,L;M - 1,M), obtained from B ljb (M;L,L) 
by removing the external unit protuberance, must be subcritical. Then, by virtue of 
Proposition 4.2, one can say that it must necessarily be (M — 1) — 2 < I* — 1, that is 
M < I* + 2. This is an absurd (see inequalities (3.17)). Then we can conclude that it does 
not exist a special saddle B\^{M; L, L) such that L + 2 > M. 

Case (ii) L + 2 < M: the internal rectangle L x L must be contained in the rectangle 
L x (m*(L) — 1), otherwise the birectangle R(L, L; M — 1, M) would be supercritical. Now 
we transform the special saddle B\^{M;L/L) into C2, a {L + 1) (notice that L > L* =>- 
L + 1 > L) and we show that the energy lowers. 

• Bi, b (M;L,L) -> R(L,L;M- l,M^AfTi = —[2 J - (h - A)]; 

• R(L,L;M -1,M) -> i?(L + 1, L; M - 1, M), A^ < because L > L > L*; 

• i?(L + 1, L; M - 1, M) -> i?(L + 1, m*(L + 1) - 1; M - 1, M), A^ 3 < because 
L > L* and L < m*(L + 1) - 1; 

• R(L + 1, m*(L + 1) - 1; M - 1, M) -> i?(L + 1, m*(L + 1) - 1; L + 2^m*(L + 1) + 2), 
AH 4 < because the external rectangle is subcritical and L + 2 < M — 1; 
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• R(L + l,m*(L + l)-l;L + 2,m*(L + l) + 2)^C 2 , a (L + l), AH 5 = 2 J - (h + A). 
Hence, 

H{B hb {M; L,L))> H(C 2 , a (L + 1)) > H{Vi) (5.10) 
for every possible choice of the dimensions M, L and L. 

Case B2. 

We consider the special saddle B 2 (M; L, L) in Fig. 5. 6. Two possible cases must be 
considered. 
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Fig. 5. 6 The internal dimensions L and L are such that the birectangle obtained by removing 
the zero unit protuberance is subcritical. The external dimensions M* — 1 and M are such that 
M>M*. 

Case (i) M* > M: the internal rectangle is subcritical, hence by removing it we obtain a 
configuration at lower energy. Then by means of arguments similar to those used in the 
case of standard Ising model (see e.g. [NS1]), one can prove that 

H(B 2 (M;L,L))>H(V 2 ) , (5.11) 
where the equality stands iff £> 2 (M; L, L) =V 2 . 

Case (ii) M* = M: see the discussion about the special saddle Bi^M; L, L). 
Case B3. 

We consider, now, the special saddle Bs(M, M; L) in Fig. 5. 7. One can easily prove that 
H(B 3 {M, M; L)) > H(d(L)) by virtue of the inequalities M < M* and L + 2 <M<M*. 
Hence, we conclude that 

H(B 3 (M, M; L)) > H(C\{L)) > H{V{) (5.12) 



B 2 (M;L,L) 
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Fig. 5. 7 The internal horizontal dimension is /*— 1 and the vertical one L is such that 
1* <L<M— 2 and L<m* (1* — 1). The external vertical dimension is M and the horizontal one 
is M<M. 



— OAAAAAAAAAAAAAAAAA 9AAAAAAAAAAAAAAAAA 

— r 1 r 1 

i — 9AA r 1 1 

I BAH 1 1 AAA — I 1 BAAAAAAA — I 

— r-r- — r — i r r — r — i 

i i' — r — i r r — r — i 

I I' — i — i r r — IIAA I 

B\ a (M, M: /.) i r — r i r r — bae r B l b ( M. M; L) 

i i' — r — i r r — r — i 

— r— r- — r — i r r — r — r 

— r— r- — r — i r r — r — i 

— r— IIAAAAAAS — r r-TIAAAAAAS- — r 

i 1 r 1 

i 1 r 1 

— nAAAAAAAAAAAAAAAAS IIAAAAAAAAAAAAAAAAS 

M M 



Fig. 5. 8 (a) The internal horizontal dimension L is such that L<L<1* — 1. The vertical one is 
L=m* (L) — 1 and it is such that L+3<M. The external vertical dimension is M and the horizontal 
one M is such that M<M. (b) The external dimensions are like those in (a). The internal 
horizontal dimension L— 1 is such that L<L— 1<1* — 2. The vertical one is L=m*(L) and it is such 
that L+2<M. We remark that for certain choices of the parameters h and A the configurations 
in (a) and (b) cannot be considered; it could be, indeed, m*(L)>M— 2. 

for every possible choice of the dimensions M, M and L. We observe that in (5.12) the 
equality holds iff B 3 (M, M; L) = V\. 

Case B4- 

We consider the special saddles B^ a {M, M; L) and B^ h [M,M\L) in Fig.5.8. First of 
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all we observe that 



H(B 4 , a (M, M; L)) > C 2 , a (L) and H(B±, b (M, M;L))> C 2 , b (L) , (5.13) 

for every possible choice of M, M and L. Equation (5.13) is a consequence of the fact that 
M < M < M* and L + 3 < M < M* in both cases. The equalities are verified in (5.13) 
iff B 4l «(M, M; L) = C 2 , a (L) or B 4 , b (M, M; L) = C 2 , b (L). 

Now, by arguments similar to those used in the discussion of case C2 we can prove that 
H(B4, a (M,M;L)) > H(Vi) and H(B 4 , b (M, M; L)) > H(Vi). 

Case B5. 

We consider the special saddle B${M; L, L) in Fig. 5. 9. 



— 9AAAAAAAAAAAAAAA 

— r r 

— r~6AAAAAAA — T 

— r-r — r— r 
„„ r „ r — r „„ r 

„„ r „ r — r — iiaa 

a -r-r — r~-eAS 

/.. /.) r |- 1 r 

--r-r- — r--r 
--r-r — r-r 
— r-nAAAAAAs — r 

— r 1 

~-r 1 

— IIAAAAAAAAAAAAAAS 

1* + 1— 

Fig. 5. 9 The internal horizontal dimension L is such that L* <L<1* — 1 and the vertical one L 
is such that M— 2<L<m* (1* — 1). The external vertical dimension is M. We remark that this 
special saddle does not exist, if we choose the parameters h and A such that m* (1* — 1)<M— 2. 



Now we transform the special saddle B$(M; L, L) into Ci(L) and we show that the 
energy lowers. 

• B 5 (M; L, L) -> i?(L, L; /* + 2, M), AH ± = -(h - A)(M - 1); 

• i?(L, L; /* + 2, M) — > - 1, L; /* + 2, M), Aif 2 < because L > M — 2 > L* and 
L < T - 1; 

• #(/* - 1, L; I* + 2, M) — > - 1, L; T + 2, L + 2), Aif 3 < since I* + 2 < M*; 

• - l,L;/* + 2,L + 2) -> Ci(L), Aif 4 = 2J-(h + A). 

By a direct calculation it can be proved that AHi + AH4 < 0, then H(B§{M; L,L)) > 
H{C\{L)) > H(V\). 
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Fig. 5. 10 (a) The internal horizontal dimension L is such that L>L*. The vertical one is 
L=m*(M— 2) — 1 and it is such that L+3>M. The external vertical dimension is M and the 
horizontal one M is such that L<M—2<1*—1. (Jo) The internal horizontal dimension L is such 
that L>L*. The vertical one is L=m*(M— 1) and it is such that L+2>M. The external vertical 
dimension is M and the external horizontal dimension M is such that L<M— 1<1* —1. We remark 
that for certain choices of the parameters h and A the configurations in (a) and (b) cannot be 
considered. 

Case B6. 

We consider the special saddles £>6, a (Af, M; L) and B§^{M, M; L) in Fig. 5. 10. 
Now, we transform the special saddle £> 6ia (M, M; L) into C 2 , a {M — 2) and show that 
the energy lowers. 

• £ 6 , a (M, M; L) -> B 6 , a (M, L + 3; L), AHx < since M < I* + 1 < M*; 

• B 6:a (M,L + 3;L) -> C 2 , a (M-2), AH 2 < because L > M-3 > I* > L* . 
Hence, we conclude that H(B 6 , a (M, M; L)) > H(C 2 , a (M - 2)) > H(Vi). 

The special saddle Be,b(M : M; L) can be transformed into the configuration C 2 ,b{M — 1) 
lowering the energy. 

• B 6)b (M,M;L) -> R(L,L;M+ 1,M), AfTi = -(/i - A)(M - 1); 

• i?(L, L; M + 1, M) -> i2(L, L; M + 1, L + 2), AH 2 < since M + 1 < Z* + 1< M*; 

• i2(L, L; M + l, L + 2) -> R(M-2, L; M + l, L+2), Aif 3 < because L > M-2 > L*; 

• R(M - 2, L; M + 1, L + 2) -> C 2 , b (M - 1), Aif 4 = 2J - (h + A). 

It is easily seen that A#i+A# 4 < 0, hence H(B 6 , a (M, M; L)) > H(C 2jb (M-l)) > H(Pi). 

This completes the proof of Lemma 5.1 in the case S < We suppose, now, 5 > ^j^- 

and observe that in this case 

m*(/*-l) =/* , (5.14) 
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as it follows from equation (5.5). In the sequel we will analyze all the cases that have to 
be discussed with arguments different from those used before. 

Case CI. 

The special saddle Ci{m) with /* < m < m*(l* — 1) — 1 cannot be considered, since 
m*(Z* - 1) - 1 = Z* - 1 (see (5.14)). 

Case C2. 

We proved above that C 2 ^ a (l* — 1) is the special saddle with lowest energy among C 2 , a (L) 
and 6*2,6 (L). This result is not dependent on the value of the real number 5. Hence, one 
can say H(C 2jb (l)) > H(C 2 , a {l)) > H(C 2 , a (l* - 1)) = H{Vi) (we remark that in the case 
8 > the special saddle C 2ya {l* — 1) and the global saddle Pi coincide). 

Case Bl. 

In order to prove that i7(£>i >a (M, M; L)) > H(Vi) we have to consider two different 
cases. 

Case (i) L > I* — 1: we transform the special saddle £>i >a (M, M; L) into V\ and we prove 
that the energy lowers. 

• B ha (M, M; L) -> R(L* - 1, L; M, M), A#i = -[2 J - (ft + A)]; 

• R(L* - 1,L;M,M) -> #(L* - 1,/* - 1;M,M), A# 2 < since L* — 1 < L* and 
L > I* - 1; 

• #(L* - 1, Z* - 1; M, M) -> #(Z* - 1, Z* - 1; M, M), Atf 3 < because Z* - 1 > L*; 

• #(Z* - 1, I* - 1; M, M) -> #(Z* - 1, Z* - 1; Z* + 2, Z* + 1), Atf 4 < since M < M* 
and M > Z* + 2; 

• #(Z* - 1, Z* - 1; Z* + 2, Z* + 1) -> Pi, A# 5 = 2J - (h + A). 
We conclude that H(B ha {M, M;L))> H(Vi) since ^)?=i A ^ < °- 

Case (ii) L* < L < I* — 1: first of all we notice that this case can be considered only if 
Z* — 1 > L*. Now we transform Bi, a (M, M; L) into Pi. 

• Bi, a (M, M; L) -> i?(L* - 1, L; M, M), A#i = -[2 J - (ft + A)]; 

• R(L* — 1, L; M, M) -> i?(Z* — 1, L; M, M), A# 2 < since L* — 1 < Z* — 1 and L > L*; 

• #(Z* - 1,L;M,M) -> #(Z* - 1,Z* - 1;M,M), A# 3 < because Z* — 1 > L* and 
L < Z* — 1; 

• #(Z* - 1, Z* - 1; M, M) -> #(Z* - 1, Z* - 1; Z* + 2, Z* + 1), Atf 4 < since M < M* 
and M > Z* + 2; 

• #(Z* - 1, Z* - 1; Z* + 2, Z* + 1) -> Pi, Atf 5 = 2J - (ft + A). 
Also in this case we conclude that H(Bi, a (M, M; L)) > H(Vi). 
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Finally, with arguments similar to those used in the case S < one can show that 
H(B 1 , b (M; L,L))> H(C 2 , a (L)) > #(7>i). 

Case B3. 

This case cannot be considered because the inequalities I* < L < m*(l* — 1) cannot be 
verified (see (5.14)). 

Case B5. 

This case cannot be considered because the inequalities M — 2 < L < m*(l* — 1) 
cannot be verified. Indeed, from (3.15) one has /* + 3 < M; hence I* < M — 2. Finally, 
m *(l* - 1) = I* =>. m *(/* - 1) < M - 2. 

The proof of Lemma 5.1 is now complete. □ 
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Section 6. The set Q and the minimum of the energy on dQ. 

In this section we define a set Q of configurations which will play a basic role in the proof 
of our results. Q will constitute an "upper estimate" of the generalized basin of attraction 
of —1, in the sense that every subcritical configuration, that is a configuration a such that 

lim P^r-i < r+i) = 1 , (6.1) 

v. 1 1 1 belong to Q; moreover given any r\ G Q there exists a downhill path leading to a 
configuration a satisfying (6.1). On the other hand there are configurations 77 G Q which 
are supercritical in the sense that 

lim P„(r_i > r +1 ) = 1 . (6.2) 

I lif crucial property of Q will be that the minimum of the energy in its boundary dQ 
will be given by V\ or Vi- 

We will see that this implies that for every configuration a with sufficiently low energy 
(H(a) < min{i/(P 1 ), J ff(P 2 )}) (6.1) is verified. 

Let us now give an example of a configuration belonging to Q which is potentially 
supercritical in the sense that (6.1) fails. 

Consider an acceptable configuration 77 which is different from — 1 only in a square 
Ao := Al (rj(x) = -1 Vx G A \ Ao); the even integer Lq will be chosen later on. 

Consider the four sublattices of spacing 2 into which Z 2 is partitioned and write 

z 2 = z\ = z 2 , a u z 2 i6 u z 2 , c u z 2 )d 

where the subscript 1 or 2 denotes the spacing; a, 6, c, d label the four sublattices of 
spacing 2. 

Suppose that Z 2 b and Z 2 d belong to the same sublattice of spacing \/2 and set 

f+1 VsGZ^nAo 

-1 Vx G z 2 6 n A 

-1 VxGZ 2 )d nA 

[ Vx G z 2 c n A 

(see Fig. 6.1). 



rj(x) 
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Starting from our definition of Q we have first to transform the — 1 with a plus spin 
among its nearest neighbours into 0. In this way we get the configuration 771 depicted in the 
left hand side of Fig.6.1. Eventually, we get a configuration fj with a unique plurirectangle 
with external edges with length L + 1 and many non-interacting unit squares of pluses 
in its interior. 

On the other hand, starting from rj we can change the — 1 in A with two plus spins 
among their nearest neighbours into +1 by decreasing the energy; we obtain the configura- 
tion r]2 depicted in the right hand side of Fig.6.1. Subsequently, still decreasing the energy 
the configuration r] 2 can be transformed into the configuration rj* = C(L — 1, Lq — 1). 

Now, if L is chosen such that /* + 1 < L < M* — 1 we have that f) (and so rf) is 
subcritical and then it belongs to Q, but rf (to which we arrived starting from 77 with a 
downhill path) is supercritical. 
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. I) . II • II • I) + - + - + - + - + - + - + -_ 

() (I I) (I I) I) I) I) I) I) I) (I 



Fig.6.1 



To construct Q, first of all we define a map T : a — > a = Ta with a an acceptable 
configuration and a a local minimum of the energy, such that the two following properties 
are satisfied 

H(a) < H(a) 

; (6.3) 

a -< a 

that is the local minimum a is bigger than a and at a lower energy level. Then we define 
the set Q as the set of configurations a such that a is subcritical, that is Pg.(r_i < r+i) — > 1 
as (3 — > 00. 

Now we define the map T : a — > <r; the definition is given in the following five steps. 
Let a be an acceptable configuration: 

{%) starting from a we construct the configuration a± by turning into zero all the minus 
spins of a which have at least one plus spin among their nearest neighbour sites. We 
remark that H{a\) < H(a) (see Fig. 3.1) and a -< o\. 
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(ii) Let us denote by c\ the minus spins cluster in the configuration <j\ which is winding 
around the torus and by c~ all the other minus spins clusters in o\. In o\ there is no direct 
interface H — , then we can conclude that every c~ cluster is inside a zero spins cluster (see 
Fig. 6. 2). Now we consider the configuration a<i obtained from <j\ by turning into zero all 
the minus spins in all the clusters c~ . The result a\ -< 02 is obvious. We have, also, that 
H{(T2) < H(cti); indeed in every cluster c~ there is at least one minus spin with two zero 
spins among its nearest neighbours; this spin can be transformed into zero lowering the 
energy. We can repeat this argument until all the spins of the starting cluster c~ have 
been transformed into zero. 
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r r r — 0AS r-nAA 

1 IAAAAA 1 1 AAA 9AE — 1 I — I 

1 llAAAE — 9 AS 1 9AE 

nAAAA r-'-'"9AAAAAS'"r 

r IIAAAAAAAS- — 0-r 

1 1 AAAAAAAAAA 9AAAAAAH 

r—0""9AAAAAS 

HAAAAAAAAAE 



Fig. 6. 2 



(Hi) In 02 there is no direct interface H — , then we observe that every cluster of plus 
spins is inside a cluster of zero spins; it can happen that in some of the plus spins clusters 
there are one or more clusters of zero spins (see Fig. 6. 3). We construct the configuration 
(T3 by removing all these clusters of zero spins. With arguments similar to those used in 
step (ii) one can prove that H(a^) < H(a2) and 02 -< 03. 

(iv) The configuration as is made of a minus spins cluster which is winding around the 
torus, the zero spins clusters denoted by Vz G {1,2, ...,/c } and the clusters with plus 
spins c+- Vz G {1, 2, k } and Vj G {1, 2, kf}. The clusters c+- Vj G {1, 2, kf} are 
all inside the cluster c°. We consider, now, the rectangular envelopes = R(c° t ) Vz G 
{1, 2, k } and the configuration cr 4 obtained by filling all these rectangles with zero 
spins; in this step the plus spins are not changed. It is immediate that H(a 4 ) < H(a 3 ) 
and cr 3 -< 04. 
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Fig.6.3 



(v) Apart from the plus spins cluster, the configuration 04 is made of zero rectangular 
clusters placed in the "sea" of minus spins. We obtain the configuration 05 by means of the 
chain construction used in [KOI], applied to the rectangular clusters R® Wi G {1, 2, k }. 

Let us briefly describe this construction. Given a set of rectangles R±, . . . , Rf we parti- 
tion it into maximal connected components cj 1 ^ with j = 1, . . . , k^ called chains of first 
generation 

(i2?...i2?) = (C{ 1) ...c2l ) ) . 

The notion of connection is given by pairwise interaction: a set . . . , R^ of rectangles 
is connected if it cannot be divided into two non-interacting parts. 

Now consider the k^ rectangles R(C^) obtained as rectangular envelope of the union 

of the rectangles belonging to Cj • Partition this set of rectangles into maximal connected 

(2) (2) 

component: in this way we construct the chains of second generation C\ , . . . We 
continue in this way up to a finite maximal order n such that the chains of the n-th 
generation are non-interacting rectangles (see [KOI] for more details). 

We call (T5 this configuration containing these non-interacting rectangular clusters 
Vz G {1, 2, /c '^} of zero spins placed in the minus spins "sea". With usual arguments 
one can prove that H(a§) < H(o~±) and o^. 

(vi) By repeating the operations described in points (iv) and (v) for the plus spins 
clusters lying in every rectangle R® Vi G {1, 2, /c '^}, we obtain the final configuration 
a. This configuration is made of the external rectangular zero spins clusters R® Vz G 
{1, 2, k '?} and the internal non-interacting plus spins clusters Vz G {1,2, k°' f } 
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and Vj G {1, 2, }. As usual one can prove that < H{a^) and (75 -< 6". 

The definition of the map T is now complete, we observe that a is a local minimum and 
that the properties (6.3) are satisfied. Finally we remark that the map T is monotone in 
the sense that 

a -< rj =>- a < fj , (6.4) 

for every couple of acceptable configurations a and r\. 
Now we state the following 

Proposition 6.1. 

U(Q) c V . (6.5) 

Namely the set of minima of the energy in the boundary of Q is contained in V . 
Proof. 

In order to prove Proposition 6.1 we consider a configuration rj G dQ and we show that 
there exists a special saddle fj such that H(rj) > H{fj). Then Proposition 6.1 will follow 
from Lemma 5.1. 

Let us consider n G Q; there exists a configuration a = n x,b with x G A and b ^ rj(x) 
such that a G Q. By virtue of the monotonicity of the map T (see equation (6.4)) and of 
the fact that a is a subcritical local minimum, it follows that b < rj(x); hence we also have 
that b ^ +1. 

We denote by R?(&) Vz G {1, k°(a)} and by #+ (<t) Vj G {1, ...k+(a)} and 
Vz G {1, /c°(cr)} the rectangles respectively of zeros and pluses which appear in the con- 
figuration a; we remark that all the rectangles R^j(a) Vj G {1, kf{o)} are inside the 
zero rectangle R®{cr). In the following, by abuse of notation, we will also denote by R®{o) 
what we will call structure R®(&), namely the complex given by the "external" rectangle 
togheter with all its "internal" rectangles of pluses (what before we called plurirectangle 
is indeed a configuration containing a unique structure). 

Case 1: b = — 1 and r](x) = 0. 

From the definition of the map T easily follows that necessarily x lies outside the 
rectangles R®{cr). 

Given the configuration f) we denote by R^ifi) Vz G {1, k°(f))} and by Rf-{fii) Vj G 
{1, kf{fj)} and Vz G {1, k°(fj)} the rectangles respectively of zeros and pluses which 
appear in it. We denote by R°(f)) the supercritical structure among the R®{fj) and by 
R\,...,R° S the rectangles of zeros such that Vz G {l,...,s} R® appears in a and Vz G 
{1, s) R\ is "inside" the rectangle R°(rj). 
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We consider, now, the configuration 771 denned as follows: rji(x) = 0, all the other spins 
are minus except for the zeros and the pluses of the structures R® Vz G {1, s}. It can 
be easily proved that H(rj) > H{r}\). We distinguish the two cases 1.1 and 1.2. 

Case 1.1: all the rectangles of pluses which appear in 771 are subcritical. 

We consider the configuration 771.1 obtained from 771 by changing into zeros all the plus 
spins. We remark that H(r)i) > if (771.1) because 771.1 has been constructed by removing 
subcritical rectangles of pluses. 

With an Ising-like argument (see e.g. [KOI]) one can prove that #(771.1) > H(J ) 2)- 
Hence in the case 1.1 we have found a special saddle with energy lower than the starting 
configuration 77. 

Case 1.2: in 77! there exists at least one supercritical rectangle of pluses. 

We consider the configuration 771.2 obtained from 771 by removing in every structure 
Vi G {l,...,s} all the subcritical rectangles of pluses and by filling with pluses the 
rectangular envelope of the union of the supercritical rectangles of pluses. We remark that 
every structure R® Vz G {l,...,s} in 771.2 is either "empty" (with no rectangle of pluses 
inside) or it has just a rectangle of pluses inside and this rectangle is supercritical. 

We denote by Q the unit square centered at the site x G A; we distinguish the two 
following cases: 

Case 1.2.1: one of the structures R® Vz G {1, s} of 771.2 (we denote it by -R1.2.1) interacts 
with Q and the structure obtained by filling with zeroes the rectangular envelope of -R1.2.1U 
Q is supercritical. 

Let us denote by 771.2.1 the configuration obtained by removing in 771 all the structures 
R® Vz G {1, s} except for R® 21 . If Q is adjacent to -R?. 2 .i then 771.2.1 is a special saddle. 
Otherwise Q is at distance one from one of the sides of the rectangle R® 2 .i or Q an d -R1.2.1 
touch in a corner; in this case it can be esily found a special saddle with energy lower than 
H(vi.2). 

Hence in the case 1.2.1 a special saddle with energy lower than the starting configuration 
77 has been found. 

Case 1.2.2: the condition 1.2.1 is not fulfilled. 

By an argument similar to the one used in [KOI] (see pages 1136-1137 therein) we can 
find two structures Ri and R2 such that: they are both subcritical, their external rectangles 
are interacting, the structure obtained by filling with zeroes their rectangular envelope is 
supercritical and #(771.2) > H(Ri) + H(R 2 ) (when we say H(Ri) with z G {1, 2} we are 
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referring to the energy of the configuration obtained by plunging the structure Ri in the 
"sea" of minuses). We still have to distinguish between two possible cases. 

Case 1.2.2.1: both structures H(Ri) with i G {1,2} have a supercritical rectangle of 
pluses inside. 

Now we consider a just-supercritical frame whose external rectangle is contained in the 
rectangular envelope of the union of the two external rectangles of R\ and of i?2- Such a 
frame surely exists and we denote it by C. 

Starting from R\ and R<i and recalling that these structures are subcritical, one can 
construct two other structures, 9ft i and 3?2 (birectangles or frames), such that the three 
following conditions are satisfied: i) #(i?i) > #(9?i) and H(R2) > #(9^); H) if the two 
external rectangles of the two structures 9ft i and 9ft 2 touch by a corner then the rectangular 
envelope of the union of the external rectangles of 9ft i and of 9ft 2 coincides exactly with the 
external rectangle of the frame C\ Hi) at least one of the two internal rectangles of pluses 
(the one in 9ft 1 or the one in 9ft 2 ) is supercritical. 

If one considers the configuration 771.2.2.1 obtained by plunging the structures in the 
"sea" of minus spins such that the external rectangles of 3ft 1 and of 9ft 2 touch by a corner, 
one can easily convince himself that #(771.2) > #(771.2.2.1)- 

Finally, starting from 771.2.2.1 we construct the special saddle 77 by performing the fol- 
lowing steps: i) we fill of zeroes the rectangular envelope of the union of the two external 
rectangles of zeroes in 771.2.2.1; H) we let grow the internal supercritcal rectangle of pluses 
until the frame C is reached; Hi) we transform into zeroes all the pluses, except for one, 
of one of the four sides of the internal rectangle, such that a special saddle is obtained. It 
can be easily proved that #(771.2.2.1) > H(fj) by comparing the energy differences involved 
in the three steps described above. We remark that the energy increase of the third step 
is largely compensated by the energy decrease involved in the second step. 

Case 1.2.2.2: one of the structures #(i?j) is "empty", in the sense that it has no rectangles 
of pluses inside. 

This case can be discussed with arguments similar to those used in the Case 1.2.2.2. 

Case 2: b = — 1 and 77 (a;) = +1. 

Starting from a one can always construct a configuration 772 such that: i) 772 G dQ; ii) 
3y G A such that 772(7/) = and 77^' 1 G Q. In this way the proof has been reduced to the 
Case 1. 

Case 3: 6 = 0. 
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The site x is inside one of the rectangles of zeroes R®(cr) Vz G {1, k°(a)}; we denote 
it by _R°. There are two possible cases that must be considered. 

Case 3.1: x is not on one of the boundary slices of _R° (the tipical situation is depicted in 
Fig.6.4). 

In this case the rectangles of zeroes in fj coincide with those in a, but the structure 
R°(fj) is supercritical (R°(fj) is the structure of fj such that its external rectangle of zeroes 
coincides with R°). 

We consider, now, the configuration 773.1 defined as follows: 1) r/3.1 is obtained starting 
from 6", by removing all the structures R\(o) Wi G {1, k°(a)} except for the one whose 
external rectangle coincides with the external rectangle of the structure R° (we denote this 
structure by R°(a)); ii) r/ 3 .i(x) = +1. It can be easily proved that H(rj) > i7(r/3.i). 

6AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

r 0AAAAAAAAA- — r 

1 BAAAAAAAAAAAAA 1 1 — I 

r ,p r nAAAAAAAAS r 

[ ]• 1 x 1 

1 1IAAAAAAAAAAAAE 0AAAAAAAAA 1 

p r p p 

p p p — 0AAA-T 

r nAAAAAAAAs — r~r~r 

r r-r-r 

r nAAs-r 

nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAS 

Fig.6.4 

We denote by R± the rectangular envelope of the union of the supercritical rectangles 
of pluses inside R (a) and by R2 the rectangular envelope of the union of the supercritical 
rectangles of pluses inside R (fj). We remark that the two sctructures R (a) and R (fj) 
have different internal rectangles of pluses, even though their external rectangles of zeroes 
coincide. 

Now we observe that there exists a rectangle R3 contained in Ri and containing R\ 
such that the configuration with all the spins minus except for the zeroes in the rectangle 
-R and the pluses in R 3 is an almost-supercritical configuration. We consider the special 
saddle fj obtained by properly putting a unit plus protuberance to one of the four sides 
of the internal rectangle of pluses of the almost-supercritical configuration found before. 
It can be easily shown that H(i]s.i) > H(fj). Hence, even in this case, we have found a 
special saddle with energy lower than the energy of the starting configuration 77 G dQ. 
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Case 3.2: x is on one of the boundary slices of _R° (see, for example, Fig. 6. 5). 

We construct the configuration 773.2 starting from a and by turning into zero only the 
spin minus at a site nearest neighbour to x. One can easily convince himself that H(fj) > 
H(V3.2)- If V3.2 £ dQ then the proof is reduced to Case 1; if 773.2 G Q the proof is reduced 
to Case 3.1. 

The proof of Proposition 6.1 is now complete. □ 



6AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

1 0AAAAAAAAA 1 

1 BAAAAAAAAAAAAA 1 1 — I 

r , r r nAAAAAAAAS r 

1 1 1 r 

1 HAAAAAAAAAAAAE HAAAAAAAAA 1 

1 1 1 1 

1 1 1 — 9AAA-T 

r nAAAAAAAAs — r-r-xr 

1 1 1 1 

r nAAs-r 

nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAS 



Fig.6.5 
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Section 7. Proof of the theorems. 

Let us first give some definitions extending the ones given in Section 4. 

We recall that by C(Zi, Z 2 ) we denote the set of configurations containing only a frame 
with internal sides Zi, / 2 . We recall the notation / := min{/i, Z 2 }, m := / 2 }. 

We denote by S(l\, / 2 ) the set of configurations obtained from C(l\, I2) by substituting 
one of the smaller internal sides with a unit square protuberance namely by substituting 
all but one plus spins adjacent from the interior to one of the internal sides of length I 
with zeroes (see Fig. 7.1). 

We denote by R(h, h) the set of configurations containing a unique birectangle obtained 
by erasing the internal unit square protuberance from >S(/i,/ 2 ) (see Fig. 7.1). We denote 
by G(/i,/ 2 ) the set of configurations obtained from the frame C(Zi,/ 2 ) by adding a unit 
square spin protuberance to one of the longer external sides in C(li, Z 2 ). A particularly 
relevant case will be the one |/i — /2I < 1 where either m = I + 1 or m = I. We remark 
that G(l — 1, /) is obtained from the birectangle R(l, I) by substituting one "free" external 
row (or column) of zeroes of length / + 2 with a unit square protuberance (see Fig.7.1); 
similarly G(l — 1, / — 1) is obtained from R(l — 1, /) by substituting one free external row 
or column of spin of length / + 1 with a unit square protuberance. 

Finally let us denote by R(h, I2) '■= -R(0, 0, l\, I2) Ui?(0, 0, Z 2 , h) the set of configurations 
without plus spins where the zero spins are precisely the ones contained inside a rectangle 
with sides equal, respectively, to Zi, Z 2 - 

We want to prove now Theorem 1. 

Let V be the set of protocritical saddles or special minimal saddles. 

If < 2A < h: V = V2 in Fig. 5.1; namely V is the set of configurations with no pluses 
and a unique cluster of zeroes given by a rectangle with sides M*,M* — 1 with a unit 
square protuberance attached to one of its longer sides. 

If < A < h < 2A and 5 < ^ then V = Vi, a ■= S(l*, I*). If < A < h < 2A and 
5 > ^ then V = V hb := S(l* - 1, /*) (see Fig.5.1). 

Now we notice that the set 5 C Oa, defined in Section 6 satisfies the following properties: 

1. g connected; -1 e Q, +1 £ Q. 

2. There exists a path u> : —1 — > V, contained in Q, with 



H(a) < H(V) Vct eu, a^V 



(7.1) 



and there exists a path uj' : V — > +1, contained in Q c , with 



H(a) < H(V) VffGw', a ± V 



(7.2) 
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O/./) S(l,l) R(l,l) 

— 6AAAAAAAAAAAAAAAAAA — 9AAAAAAAAAAAAAAAAAA — 6AAAAAAAAAAAAAAAAAA 
-T-9AAAAAAAAAAAAAAT— T-9AAAAAAAAAAAA— r — r-9AAAAAAAAAAAA— r 

1 i- 1 r r-r 1 r r-r r— r 

1 i- 1 r r-r iiaa r i r r— r 

1 i- r-r— -r-r 9As-r n r— r 

1 i- r-r— -r-r r -r -r-r r-r 

r , r r-r -r-r r -r -r-r r-r 

1 i- r-r -r-r r-r-— r-r r-r 

r-nAAAAAAAAAAAAAS-r— -r-nAAAAAAAAAAAS— r— -r-nAAAAAAAAAAAS— r 

— IIAAAAAAAAAAAAAAAAAS — IIAAAAAAAAAAAAAAAAAS — IIAAAAAAAAAAAAAAAAAS 
/ + 2 / + 2 / + 2 
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— 6AAAAAAAAAAAAAAAA- — 9AAAAAAAAAAAAAAAA 6AAAAAAAAAAAAAAAA 

-T-9AAAAAAAAAAAAT r-OAAAAAAAAAAAAT r-9AAAAAAAAAAAAT 

1 I- 1 r — r-r r-r — r-r r-r 

1 i r-r- — r-r r-r — r-r r-r 

r , r 1 iiaa r-r r-r — r-r r-r 

1 i r-9AH— r-r r-r — r-r r-r 

r , r 1 |- — r-r r-r — r-r r-r 

1 i- r-r- — r-r r-r — r iiaaaa 9AAAAAS r 

r-nAAAAAAAAAAAST r-IIAAAAAAAAAAAST r — IIAS T 

IIAAAAAAAAAAAAAAA5- — IIAAAAAAAAAAAAAAAS IIAAAAAAAAAAAAAAA5 

/ + i / + i / + i 



/ t > ( / i. /) G(l-1, l-l) C(l-1, 1 1) 

9AAAAAAAAAAAAAAAA- — 9AAAAAAAAAAAAAAAA- — 6AAAAAAAAAAAAAAAA— - 

•— r-6AAAAAAAAAAAAT r-0AAAAAAAAAAAAT r-9AAAAAAAAAAAAT 

1 i- 1 |- — r-r 1 r — r-r r-r 

1 i- 1 |- — r-r r-r — r-r r-r 

r , r 1 |- — r-r r-r — r-r r-r 

r , r 1 |- — r-r r-r — r-r r-r 

1 i- r-r- — r-r r-r — r-r r-r 

r-nAAAAAAAAAAAST r-IIAAAAAAAAAAAST r-IIAAAAAAAAAAAST 

1 r HAAAAAAA 0AAAAAAS IIAAAAAAAAAAAAAAAS 

IIAAAAAAAAAAAAAAAS HAS 

/ + i / + i / + i 



Fig. 7. 2 Contraction of a squared frame. The energy differences in- 
volved in each single step of the contraction are: (h+\)(l-i), -[2j-(h+\)], 

(h-X)(l+l), -[2J-(h-X)], (h+\)(l-2), ~[2J-(h+X)], (h-X)l, ~[2J-(h-X)]. 

In the case V = V\ (7.2) easily follows from the arguments of proof of Proposition 4.1: 
uj is constructed following a sequence of shrinking subcritical droplets whereas u>' follows 
a sequence of growing supercritical droplets. In the case V = V2 (7.2) follows from the 
arguments of proof of Proposition 4.2. 
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3. The minimal energy in dQ is attained only for "protocritical" (global saddle) configu- 
rations a &V; namely, 

mm (H(a) - ff(-l)) = H(V) - ff(-l) =: F , (7.3) 

mm (H(a) - H{V)) > . (7.4) 

aedg\v 



We notice that, starting from any a G V, we can change a spin adjacent to the unit 
square protuberance always present in V (from — 1 to in V2 if h > 2 A and from to 
+1 in V\, a or Vi t b if h < 2A) in order to get a "stable protuberance of length 2". This 
protuberance is called stable since its growth takes place decreasing the energy while its 
shrinking requires an increase of energy. The probability of the above described single spin 
change is not smaller than j^y (see, for instance, [NS1] for more details on this point). 

In other words, with probability separated from zero, uniformly in starting from V, 
we reach the strict basin of attraction of a supercritical minimum. Then, for any e > 0, it 
follows from Proposition 4.1 that the probability to reach +1 before reaching —1, can be 
bounded from below as: 

P v (t +1 < r_i) > exp(-e/3) . (7.5) 

We get from Proposition 4.1 that, for (3 sufficiently large, the typical time, starting from 
V\ to reach +1 is much shorter than the typical time to get to V starting from —1 

lim P Vl (r +1 < exp(r 1 ) | r +1 < r_i) = 1 . (7.6) 

:< 'i a suitable Ti < V. 

Moreover by an analysis totally analogous to the one needed for the Ising model (see 
for instance [NS1]) one can get the same results starting from V2] namely 

lim P V2 (t +1 < exp(r 2 ) I t+i < r_i) = 1 . (7.7) 

for a suitable T 2 < T. 

In Appendix A we state and prove a result concerning the sequence of passages through 
V and the typical time to see an "efficient" passage through V namely one followed by a 
descent to +1. 

From Propositions 3.4, 3.7 in [OS1], Proposition A.l of Appendix A, (7.5), (7.6) we 
easily get Theorem 1. □ 
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We want now to give the definition of the tube T of trajectories appearing in the state- 
ment of Theorem 2 below. It represents the typical mechanism of escape from metastability 
in the sense that, with probability tending to 1 as (3 tends to infinity, during its first ex- 
cursion from —1 to +1, our process will follow a path in T. 

T will be optimal in the sense that it cannot be really reduced without loosing in 
probability. 

T involves a sequence of "droplets" with suitable geometric shapes and suitable "resis- 
tance times" in some "permanence sets" of configurations related to these droplets. The 
precise statement about the typical paths during the first excursion between —1 and +1 
will involve a certain randomness of these resistance times inside the different permanence 
sets appearing in T. 

In T we will distinguish two parts. The "up" part T u namely the ascent to V and the 
"down" part Ta from V to +1. This second part is almost downhill in the sense that 
all the paths ui = o~o, ui, . . . , Uj, . . . e Td will be such that: 

a = V, 3 f : a T = +1, max H(a) < H(V), mm H (a) = H (+1) . 

Whereas T u is almost uphill in the sense that all the paths u = ctq, a±, . . . , . . . e T u will 
be such that: 

(7 = -I, 3 f : a? = V, max H(a) < H(V), min H(a) = H(-l) . 

aeLo\v <?euj 

In the following we give the definition of the time reversed tube T of T u . 

T will also be almost downhill; it will describe the typical first "descent" from the 
protocritical saddle to —1. By general arguments based on reversibility (see ref. [SI]), 
we will deduce the desired results on the first excursion from — 1 to V saying that with 
probability tending to one as (3 tends to oo it takes place in the tube T u . Then to conclude 
our construction of T we will only have to determine T^. 

Let us now recall some basic definitions of [OS1] concerning the first descent from any 
configuration r/o contained in a given cycle A to the bottom F(A) valid not only for our 
Blume-Capel Metropolis dynamics but also for a general "low temperature" Markov chain 
satisfying Hypothesis M in Appendix A. We refer to Appendix A where this more general 
set-up is introduced. 

We will first define in general the set of "standard cascades" emerging from a configura- 
tion rjo; our intention is to apply a (simplified version of a) result of [OS1] telling that with 
high probability when (3 — > oo the first descent from r/o to F(A) follows, in a well specified 
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way, a standard cascade. Thus the main model-dependent work will be to determine, in 
our specific case, the set of standard cascades. In particular we will reduce the problem of 
the determination of the tube of typical trajectories followed by our process during its first 
descent to — 1 starting from a configuration <ro in Q immediately reached starting from the 
global saddle V, along a downhill path entering Q, to find the set (denoted by T) of all 
the standard cascades (in a suitable cycle) emerging from do- 

A standard cascade emerging from a state r/o is a sequence: 

T(rj , ui, 7/1, u>2, . . . , tjm-1, wm) = wiUQiUw 2 ,...,Qm-iUwmUQm, (7.8) 

where for i = 1, . . . , M: Ui is a downhill path emerging from r\i-\ and ending inside the 
"permanence set" Qi. Each path u>i can be downhill continued up to a stable equilibrium 
point £j e Qi. The Q^s are special sets being a sort of generalized cycles containing also 
the minimal saddles between £j and F(A); for % = 1, . . . , M — 1 ?7j e are minimal saddles 
between & and -F(A); finally £m C Qm C -F(^) ( see [OS1], Section 4 for more details). 

Notice that u>i can just reduce to one downhill step from r/j_i to Q^; in this case we use 
the convention: Ui = rji-i- 

We do not give here the precise definition of the Q^'s since it happens that we do not 
really need it. In our particular case of Metropolis dynamics for the Blume-Capel model 
with particular initial conditions (of interest for our applications) as we will check we have 
some semplifications w.r.t the general case. 

The most important is that the Qi's, for % = 1, . . . , M — 1 are replaced by genuine cycles 
Af, rji, not contained in Ai, is an element of S(Ai) and u>i ends in the interior of Ai. 

We will apply the general theory developed in [OS1] to two cases. In the first one, 
when analyzing T the cycle A will be the maximal connected set A in Qa containing — 1 
with energy less than H(V). It follows from Proposition 3.4 in [OS1] that A is contained 
in the set Q introduced in Section 6 and that S(A) = V. Always in this case we have: 
F(A) = Q m =^m = -1. 

In the second case, in the study of 7^ the cycle A will be the maximal connected 
component A in Oa containing +1 with energy less than H(V). It is immediate to see 
that A C Q c . In this case we have: F(A) = Qm = £m = +!• 

In both cases, as we said before, for suitable initial conditions we will verify that the 
QiS for % = 1, . . . , M — 1 are replaced by genuine cycles A t ; M will depend on the initial 
configuration r/o as well as on the particular choice of the parameters J, h, A. Ui ends 
in the interior of Af, rji E <S(£i, F(A i+ i)), not contained in Ai as we said before are 
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minimal saddles in the boundary dAi. The cycles Ai are precisely the maximal connected 
components containing £j with energy less than H(rji) (£j G A± are the minima towards 
which Ui can be downhill continued). 

We consider an initial configuration t]q corresponding to one of the following five cases: 

1. A = A: rjo G A D [R(l, I) U R(l, I + 1)] for some / > L. 

2. A = A: 7] G R(M*, M* — 1). 

3. 4 = i, < A < /i < 2A and 8 > r] G B(C(l* - 1, /*)). 

4. A = i, < A < /i< 2A and 6 < r] G B(C(l*, I*)). 

5. A = A: r] e B(R(M*,M*)). 

Remark. 

We could even consider much more general initial configurations 770- It is not true (see 
the definition of the set D in Section 4 that for any r]o the simplified version (involving 
genuine cycles Ai in place of the sets Qj) of the general [OSl] results holds true. In fact 
with the very particular choice t]q G A n [R(l, I) U R(l, I + 1)] as we will see an even simpler 
statement holds: the Ui will be almost all coinciding with rji (in the above specified sense). 

Warning. 

We want to warn the reader of the use that we are going to make, in the construc- 
tion of the tube T, of the equivalence class of configurations as it has been specified in 
the remark in the proof of Proposition 4.1. In fact, strictly speaking what we will con- 
struct and call standard cascades are sets of standard cascades obtained from equivalence 
classes of configurations modulo rotations, translations, inversion and "displacement of 
protuberances" . 

Let us now start with the definition of the set T of the standard cascades emerging 
from a configuration ctq in Q immediately reached starting from the global saddle V, along 
a downhill path entering into Q. 

We consider first the case a = h/X < 2 . The other case of a > 2 is almost identical to 
the corresponding one for the Ising model and will be treated later on. 

We have to distinguish two cases: 5 = I* — 2J ~^ t ~ A - > < A^p, when the global saddle 
V has the form V\, a = S(l*, /*) given in Fig.5.1; or 6 > when the global saddle has 
the form Vi,b = S(l* - 1, /*) also given in Fig.5.1. 

Let us first consider the case 5 < (like in Fig.7.2 for / = /*). Let V\ = R(l*,l*) 

be the configuration obtained from V\ by erasing the unit square protuberance. V\ is a 
subcritical birectangle; it belongs to the set Q and satisfies condition 1 above. 
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To construct the tube T we have basically to solve the above described sequence of 
minimax problems. To simplify the exposition we divide the tube T into four segments 
corresponding to four different mechanisms of contraction; we write: 



T = 7i U T 2 U T 3 U T 4 . 



(7.9) 



The most relevant ones are the first and the second part. As we will see the third part for 
h < 2 A reduces just to a simple downhill path. 

We start from the determination of the minimal saddle r/i := S(Vi,—l) between V\ 
and —1. 

From the results of Section 4 we know that S{V\, — 1) is not trivial in the sense that it 
differs from V\ and 



Thus the first "permanence set" Qi of our standard cascade is the cycle Ai*_ij* defined 
as the maximal connected set of configurations containing R(l*J*) with energy less than 
H(S(l* — 1, /*)). We recall that the basic inequality to be checked in order to get (7.10) is 



which is verified for L* < I < I* — 1. 

For any /: L* </</* — 1 we define the cycle A^i [Aij +1 ) as the maximal connected 
set of configurations containing R{l,l) (R(IJ + 1)) with energy less than H(S(l — 1,/)) 
(H(S(IJ))) (see Fig. 7. 2). By extending the previous argument we get that the first part 
of our standard cascade is given by: 



Vl := S{V U -1) =S(l*-l,l*) . 



(7.10) 



H(S(l,l))-H(S(l-l,l)) > 



Tx = S(l* — 1, /*), S(l* — 1, /* — 1), Ai»_ 2) Z*_l, 



...,S(L,L+1),A Z ~,S(L,L) 



(7.11) 



Then we observe that for I < L — 1, we have 



H(S(l,l)) < H(G(l,l)) 
H(S(l,l+l)) < H(G(l,l + l)) ' 



(7.12) 



(7.12) are the basic inequalities to get that, for l < I < L: 



S(R(l,l + l),-l) = G(l,l) 
S{R{l,l),-l) = G(l-l,l) 



(7.13) 
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It is clear from the results of Section 4 that the subsequent permanence sets are the cycles: 

A z,Z-i' a z,l-i' a Z-i,l-v a Z-i,l-v • • • ' A h^ A h> A h-^ A h-^ • • • > A t,i ( 7 - 14 ) 

where lo = [j + 1], /o < we notice that for our present choice of the parameters: 
A < h < 2 A, we have Z = 2 but we could consider a general situation / > 2 as well 
when analizing the contraction of a subcritical frame in the region h > 2 A (case 1 above). 
Moreover, for Iq < I < L: 

A\ l _ 1 = maximal connected component of the set of configurations containing R(l,l) 
with energy less than H(G(l — 1, /)); namely Aj t _ x is the strict basin of attraction of R(l, I): 

Al_ x = B(R(l,l)) , 

with bottom 

and minimal saddle 

=G(l-l,l) ; 

Aj t _ t = maximal connected component containing C(l — 1,1) with energy less than 
H(S(l - 1,1)). We have: 

Al_ ± = B(C(l-l,l)) 
F(Ali-i) =C(l-l,l) 

and 

SiA^) =S(l-l,l) . 

For lo + 1 < I < L we define A\_ x { _ x = maximal connected component containing 
R(l — 1, I) with energy less than H(G(l — 1,1 — 1)). We have 

=B(R(l-l,l)) 
F(A}-i : i-i) =R(l-l,l) 

and 

SiAl,^) =G(l- 1,1-1) . 

Af_ 1 1_ 1 = maximal connected component containing C(l — 1,1 — 1) with energy less than 
H(S(l- 1,1-1)). We have: 

AU,i-i =B(C(l- 1,1-1)) 
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F(Ali,i-i) =C(l- 1,1-1) 

and 

S(i4?_i,«_i) =5(1-1,1-1) • 
Then the second segment of the standard cascade is: 

,G(L-l,L),Al~ ,S(L-1,L),A± ~ , G(L — 1,L — 1), A~ ~ , 

L,L-l' v ' y ' L,L-l' v ' ; ' L-l,L-l' v ' y ' L-l,L-l' 

5(L-l,L-l),...,5(Zo,/o),i4U-i • ( 7 - 15 ) 

We notice that both the first and the second part Ti , T 2 of the tube T describe a contraction 
following squared or almost squared frames; but whereas in the first part the permanence 
sets are cycles with many minima in their interior, in the second part they are "one well" 
in the sense that they coincide with the strict basin of attraction of their bottoms. The 
typical times spent inside these cycles and the typical states visited before leaving them 
are different in the two cases of T{ and T 2 . 

The third part T3, that we are going to define, corresponds to the shrinking of the 
interior rectangle of the frame C(/ , lo)- Indeed it follows from Section 4 (see (4.6) therein) 
that for / < l the lowest minimal saddles in the boundary of the basin of attraction 

A i ,i -i of the birectan g le R ( l o, lo) = R(lo, k -lJo + 2, l + 2) U R(l - 1, / , k + 2, l + 2) 
is not G(l — l,/o) corresponding to £2 in Fig. 4. 6 but, rather, the saddle £3 in Fig. 4. 6; 
in other words starting from the birectangle _R(/o,/o) it is no more convenient to continue 
the contraction along frame shapes but, on the contrary, the internal rectangle starts its 
independent shrinking keeping fixed the external one. It appears clear that if h < 2A 
then the shrinking and disappearing of the internal two by two rectangle is just a down 
hill path where the number of internal plus spins decreases monotonically to zero. If we 
were considering a general initial condition corresponding to the above case 1 namely the 
contraction of a subcritical frame for h > 2A, then we would have had Iq > 2 and the 
shrinking and disappearence of the internal rectangle would have followed a sequence of 
squared or almost squared rectangular shape exactly like in the case of the standard Ising 
model. 

In the following we will consider birectangles R(Li, L2; Mi, M2) (see (3.3)) also for 
L U L 2 = 0,1. 

Then the third part for h < 2A is just the downhill path: 

% = SI,R(1, 2; 4, 4), R(l, 1; 4, 4)^(0, 0; 4, 4) (7.16) 
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where by S3 we denote the saddle depicted in Fig. 4. 6 when the external rectangle is a 4 x 4 
square and the internal cluster is a "triangle made by 3 sites" . 

Finally the fourth part is just an Ising-like contraction of the remaining 4x4 rectangle 
of zeroes. We will observe first a sequence of permanece sets (corresponding to stable 
rectangles) and saddles and finally the downhill path describing the disappearence of the 
last 2x2 stable rectangle. 

We have: 

T4 = S±,Ri, S 2 , R2, S3, R3, S4, _R 4 , u (7.17) 

where R x = R(0, 0; 4, 3) U R(0, 0; 3, 4), R 2 = R(0, 0; 3, 3), R 3 = R(0, 0; 3, 2) U #(0, 0; 2, 3), 
i?4 = -R(0, 0; 2, 2); the downhill path ui is given by 

Co :=S 5 ,R 5 ,R 6 ,-1 , (7.18) 

with R 5 = R(0, 0; 1, 2) U #(0, 0; 2, 1), R 6 = R(0, 0; 1, 1); the saddles S it i = 1, . . . , 5 are 
obtained from the rectangles in Ri by adding a unit square protuberance to one of its 
longer sides. 

This concludes the definition of T for h < 2 A, 5 < 4^-- 

In the case h < 2A, 5 > 4^ the definition of T is almost identical; we only have 
to modify a little bit at the very beginning the definition of Ti by eliminating its first 
permanence set. 

Indeed we know from Section 5 that now H(S(l* — 1,/*)) > H(S(l*J*)) so that the 
protocritical saddle is, in this case, = S(l* — 1,1*). Now the configuration V\ obtained 
from V\ by erasing the unit square protuberance is the subcritical rectangle V\ = R(l* — 
1,/*); again this belongs to case 1. Then the first permanence set is now Ai* and 
we have 

7i = A,..!,,..!, S(l* -1,1*- 1), A*-2,/*-i • • • S(L, L + 1),A Z ~, S(L, L) . (7.19) 

The other segments of the tube %, i = 2, 3, 4 are defined exactly as before. 

The last case that we have still to analyze to construct T is h > 2A. In this case 
the protocritical saddle is V = V2 and the tube T is just an Ising-like contraction along 
squared or almost squared rectangular clusters of zeroes in a sea of minuses. 

Now the configurations obtained by erasing the unit square protuberance, containing a 
unique subcritical rectangle of zeroes in a sea of minuses is given by: 

P 2 = R(M* - 1,M*) ; 
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notice that V2 is included in the case 2 above. 

We observe that the appearence of a single plus spin will induce the overcoming of an 
energy barrier greater or equal to 4 J — (h + A). It is very easy to see that we can proceed in 
the construction of the set of the standard cascades emerging from R(M* — 1, M*) without 
be forced to overcome a barrier larger then 2J so that certainly in all these standard 
cascades, for our choice of the parameters, we will never see a single plus spin appearing. 
Indeed one easily convince himself that the sequence of minimax problems to be solved 
are the exact analogue of the ones arising in the analysis of a subcritical contraction for a 
standard Ising model. We refer to [SI], [KOI] for more details. For completeness in the 
following we summarize the results using our notation. 

The first permanence set is B(R(M* - 1, M*)). 

Let us define the following sequences of couples of integers: 

(Zi, mi), (h, m 2 ), . . . , (l N , m N ) N = 2M* - 2 
(Zi,mi) = (M* - 1,M*), (l N ,m N ) = (1, 1); \k - rrii\ < 1 : m^I; or m i = l i J rl 

if (k, mi) = (1,1+ 1) then (l i+1 ,m i+1 ) = (/, /) 

if (k, mi) = (1,1) then (l i+1 ,m i+1 ) = (/ - 1,/) . 

Given (/, m) as before: for |Z — m\ < 1, 1 < m < M* — 1 we denote by S(l, m) the saddle 
obtained from R(l, m) by adding a unit square protuberance (with a zero spin inside) to 
one of its longest sides. 
We have: 

T = B(R(h, mi)), S(h, m 2 ),B(R(l 2 , m 2 )), . . . , S(l N , m N ), R(l N , m N ), -1 . 

This concludes the definition of T. 

Let us now pass to the definition of the descent part Td of the tube T. 
We start from the case h < 2A, 5 > . 

It is immediately seen that by adding to V\ = Vi,b = S(l* — 1,1*) a unit square 
protuberance to form a stable protuberance of length 2 we get a configuration t]o included 
in case 3. 

We distinguish in 7^ two parts: 7^,1 and Td : 2- 
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For I* — 1 < I < M — 2 we denote by Ai_ij the cycle given by the maximal connected set 
of configurations containing C(l — 1,/) with energy less than H(S(l,l)). We easily verify 
that F(Ai-ij) = C(l - 1, /), SiAi-u) = S(l - 1, /). 

For /* < / < M — 2 we denote by Aij the cycle given by the maximal connected set of 
configurations containing C(7, /) with energy less than H(S(l 7 1 + 1)). We easily get that 
F{A ul ) =C(l,l), S(A M ) = S(l,l + 1). 

For /* — 1 < I < M — 2 we denote by f^_i,j the set of downhill paths starting from 
S(l — 1, /) and ending in Ai-ij. Similarly, for /* < / < M — 2 we denote by Q^i the set of 
downhill paths starting from S(l,l) and ending in A\^. We set 

r d)1 = 7/o,Ai.-i,i., 5(r,r),n /%/ *,A*^,5(r,r + i),n^,^ +1 ,...,5(M-2,M-2) . 

As it has been shown in Section 4 for / > M — 2 the growth is typically symmetric in the 
sense that the probability of growth in the directions parallel or orthogonal to the shortest 
side of our supercritical frame are logarithmically equivalent for large (3. Moreover it 
follows from the analysis developed in Section 4 that for I > M — 2 the set T> defined in 
(4.14) do not play any particular role and the permanence sets are cycles given by the 
strict basins of attraction of frames C(h,l2) or birectangles R(h,l2) • The second part 
T dj2 of T d will describe the supercritical growth starting from I = M — 2. To construct T dy2 
we need some more geometrical defininitions. 

For a given frame C(/i, Z2), we use the notation C(Z, m) to make explicit the shorter and 
longer sides / and m, respectively. 

We denote by G>(/,m), G<(/,m), respectively, the saddle configurations containing a 
unique droplets obtained by attaching a unit square protuberance (with a zero spin inside) 
to a longer or shorter external side of C(Z, m). 

We denote by _R> (/, m) , _R< (/, m) , respectively, the birectangles obtained from G > (/, m) , 
G<(/,m) by extending the unit square protuberance to an entire side. 

We denote by S > (l,m), S < (l,m), respectively, the saddle configurations containing a 
unique droplet obtained from i?>(/,m), i?<(/,m) by attaching a unit square protuberance 
(with a plus spin inside) to the internal free side. 

We denote by fi> (/, m), fi< (/, m), respectively, the set of all the downhill paths emerging 
from Sy^^m), S < (l,m) and ending in B(C(l + 1, m)), S(C(/,m + l)); finally we denote 
by 0>(/,m), 0<(/,m) the set of all downhill paths emerging from G>(/,m), G<(/,m) and 
ending in £?(!?>(/, m)), B{R < {l,m)). 

Given (/,m), we denote by r>(/,m) the sequence: S(C(/,m)), G>(/,m), 0>(/,m), 
S(i?>(/, m)), S>(Z, m), 0>(Z, m), B(C(l, m + 1)). Similarly we denote by T<(/, m) the se- 
quence: B(C(l,m)), G<(/,m), 0<(/,m), S( J R < (/,m)), S<(/,m), 0<(/,m), S(C(/ + l,m)). 
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A sequence (ii, ?7ii)i=i,2,... with < rrii is called regularly increasing if: 

l\ = mi = M — 2 and, for any i = 1, 2, . . ., either (Zj+i, m, + i) = (Z^, := (/^ + 1, mj) 
or (Z i+1 , m i+i ) = (k, m i ) < := m s + 1). 

If li = rrii = L = the side of our torus A, we set Zj + i = m l+ i = L. 

Let £ be the set of all regularly increasing sequences. For any (hjirii)^^,... G £ we 
define: 5(l h mi) := > if (Z i+ i,m i+ i) = (/ i ,m i ) > and 5{li,mi) := < if (Z i+ i,m i+ i) = 
(l i ,m i ) < . 

From the arguments developed in Section 4 it easily follows that the second part of T d 
is given by: 

Td,2 = ^{h,m i ) i=ia ,...eC^5(h,m 1 ){h, ™l) U ^S{l 2 ,m 2 ) (^2, "7, 2 ) U . . . , 1^ . m .) (Zj, m,), ... . 

This concludes the construction of Td for the case h < 2A, 5 > 4^ . 

The case ft < 2A, S < requires only minor changes: the only difference is that 

now we have to start a step further. Indeed it is immediately seen that by adding to 
V\ = Vi :a = S(l*,l*) a unit square protuberance to form a stable protuberance of length 
2 we get a configuration 770 included in case 4. We have 

T dA = 770, A l%l * , S(l*, V + 1), • • • , >5(M* — 2, M* — 2) . 

The rest is identical. 

For the case h > 2A we have exactly the same behaviour as in the Ising model namely we 
pass to consider an initial condition like in the case 5. Then we have a symmetric growth 
along a sequence of supercritical growing rectangles of zeroes in a sea of minuses up to 
the configuration 0. Subsequently we have again an Ising-like nucleation of a protocritical 
droplet of pluses in the sea of zeroes (an L* x (L* — 1) rectangle with a unit square 
protuberance attached to one of its longer sides) up to the configuration +1. This last 
case has been already analyzed in detail (see, for instance [NS1], [SI]). We leave the details 
to the reader. 

One can easily convince himself that this indeed concludes the construction of the set of 
all standard cascades emerging from any of the above specified five type of initial conditions 
for any value of the parameters (not only for the subcases that we have explicitely treated) . 

We can now state our main result on the tube of typical trajectories during the first 
excursion between —1 and +1. 

Let T := T u U V U Td with T u given by the time reversal of the set of standard cascades 
in A C Q emerging from V: T u := 1ZT (the time reversal operator acts on paths in this 
way: for ui = (x\, x%, ■ ■ .xn-i, %n) '■ 7^ = xn-i, ■ ■ ■ x i))', T^d given by the set of 
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standard cascades in A C Q c emerging from V. Let V\ be either V\ or Vi according to 
the values of the parameters J, h, fx; let V\ be either V\ or P 2 according to the values of 
the parameters J, h, [i. 

Theorem 2. 

Consider the dynamical Blume-Capel model described by the Markov chain with transition 
probabilities given in (2.6) of Section 2. For any choice of the parameters J, h, \i compatible 
with (3.17). 



lim P_i((7 t eTVte [f_i,r+i]) = l . 

I lif history of the process in T is described in the following way: 
consider an initial configuration rjo corresponding to one of the following five cases: 

1. A = A: i] G A n [R(l, I) U R(l, I + 1)] for some / > L. 

2. A = A: 7] G R(M\ M* — 1). 

3. A = A, < A < h < 2A and 5 > r/ G B{C{1* - 1, /*)). 

4. ^ = i, < A < h < 2A and S < %±A; 770 G S(C(/*, /*)). 

5. A = A: 770 G B(R(M*,M*)). 

Then, considering for any such A, ?]o the set of all standard cascades emerging from r]o 
and falling into F(A) we have 

ii) 

3 5 > such that lim P m (r F(A) < exp(P[H(rn) - H(F(A)) - 5]) = 1 , 

,3^00 

zzz) 

lim (V t < t F ( A \ : x t G T(-q , u u i] U u 2 , ■ ■ ■ , Vm-i, wm) 

p— >oo 

for some standard cascade 770, o>i, 771 , a;2, . . . , r\M-\i ^m) = 1 , 

iv) moreover, with probability — > 1 as (3 — > 00, there exists a sequence r}o,u>i,r}i,u)2, 
. . . , rjM-i, <*>m such that our process starting at t = from 770 , between t = and 
* = r F(A) 5 after having followed the initial downhill path u)\, visits, sequentially, the 
sets Ai, A2, . . . , Am-i exiting from Aj through rjj and then following the path uJj+i before 
entering Aj+±. 

For every e > with probability tending to one as (3 — > 00 the process spends inside each 
AjZtimeTj^-.eMPiHiv^-HiFiAjV-e}) < T 3 {e) < exp(/?[7J(r ?J )-if(F(A J ))+ £ ]) 
and before exiting from Aj it visits each point in Aj at least exp f3e times . 
Proof. 



l//ebi)raio/2008 [83] 



7:83 



We easily get that 



lim P Vl (ai e Vx\o x e Q) = 1 . (7.20) 



Indeed (7.20) follows from the fact that there is only one first possible step in any downhill 
path from V\ to Q: it corresponds to erasing the unit square protuberance to get V\. 
On the other hand we have: 

lim P Plb (a 1 = V o\a 1 eg c ) = l . (7.21) 

The proof is an immediate consequence of Theorem 1, (7.20), (7.21), Theorem 1 in [OS1], 
and the results of [SI]. □ 
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Section 8. Conclusions. 

We have described the metastable behavior of a dynamical Blume-Capel model. Our 
updating rule is given by the classical Metropolis algorithm but it is clear that our results 
extend to a wide class of single-spin-flip reversible dynamics. 

Our results refer to the asymptotic regime of small but fixed magnetic field h and 
chemical potential A, large but fixed volume A and very large inverse temperature (3. 

We take mainly the point of view of the so called pathwise approach to metastabil- 
ity aiming to describe the typical behaviour of the random trajectories of our stochastic 
dynamics rather than describing the evolution of the averages. 

Blume-Capel model exhibits the interesting feature of the presence of three possible 
phases. The equilibrium phase diagram is, consequently, very reach and interesting. The 
most important aspect from the point of view of the study of the dynamics of metastability 
is the presence, near the triple point, of two competing metastable phases. This means 
that, for instance, if one wants to describe the decay from the metastable — 1 phase to the 
stable +1 phase one has to take into account the presence of another metastable phase: 0. 

We took as initial condition the state —1 and we analyzed the region of parameters 
< A < h. Let us subdivide it into the regions II and III defined as follows: 

II := < A < h< 2A 
III := < 2A < h 

It is easily seen that, with the same arguments developed in Sections 3,4,5,6 we could 
analyze the region 

IV := < -A < h 

as well. In II, III, IV the stable equilibrium phase (absolute minimum for the energy) is 
+1 and we have: 

H{-1) > H(0) > H{+1) . 

In the region 

I : = < h < A 

we have 

H(0) > H(-l) > H(+l) 

and then it is reasonable to expect and not difficult to prove that in the decay from — 1 
to +1 the state does not play any role. Indeed it is sufficient to exhibit a mechanism of 
transition from —1 to +1 involving an energy barrier smaller than H(0) — H(—l). 
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This is very easy to achieve if the volume A is sufficiently large. 

In this paper we analyzed in detail the regions II and III which happen to be, in a sense, 
the most interesting ones. In the region IV one has the same local minima for the energy 
as in the regions II and III; they are sets of non-interacting plurirectangles; but now the 
comparison between the times ^1,^2,^3,^4 introduced in (3.13) changes totally. The main 
difference w.r.t. the regions II, III is that now, in IV, we have: 

M* < L* , 

and so we cannot even consider a possible mechanism of nucleation along a sequence of 
frames. Indeed one has that a birectangle is supercritical if and only if the minimal external 
side is not smaller than M*. Then, like in region III but in a much easier way, we can prove 
that the escape from — 1 starts with an Ising-like nucleation of a protocritical droplet V2 
leading to 0. But now, contrary to the region III the typical time T~-~^- for going from 
— 1 to is much shorter than the typical time T-~^ + - for going from to +1 so that the 
asymptotics of the time T~-~^ + - of the transition from — 1 to +1 is dominated by T-~^ + -. 

The situation in which a priori one could expect a competition between the two 
metastable phases would be at a first glance the union of the regions II, III, IV. By 
arguing more carefully with a heuristic analysis of the heights of the possible barriers be- 
tween —1 and and between —1 and +1 (given by the energy of formation of suitable 
critical droplets) one is led to expect that the two metastable phases corresponding to — 1 
and are in a sense really competing only around the half line < h = 2A separating the 
regions II and III. This value h = 2A depends on the particular form of the Blume-Capel 
hamiltonian. 

The main result of the present paper consists in the rigorous proof of the above heuris- 
tics. 

From mathematical point of view we had to solve some large deviation problems. This 
kind of problems would be extremely hard for a general non-reversible dynamics but their 
treatment is very much simplified by the reversibility property of the dynamics. 

In particular to get the result we had to solve the minimax problem of the determination 
of the global saddle between —1 and +1. This is the really hard point of the work. We 
could handle the large deviation problems a la Freidlin-Wentzell arising in the study of 
some rare events in the framework of our low temperature Metropolis dynamics by taking 
advantage of a general approach to the study of typical trajectories, during the first exit 
from a non-completely attracted domain, recently developed in [OS1]. Nevertheless we 
still had to face the crucial model-dependent part consisting in solving some geometrically 
quite involved variational problems. 
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In particular we had to exclude, as highly depressed in probability, any mechanism of 
transition based on coalescence and we had to single out, among many others, only very 
few possible mechanisms of nucleation. 




h/2 



Fig.8.1 



We were able to rigorously compute the lifetime of the metastable state, namely the 
tipical transition times T\,/j, for different values of the parameters A, h. It turns out that 
these transition times are given by 

T X)h ~ exp(/?r Aj/l ) 

where the "activation energy" , for very small values of A, h has the following expression 



\,h 



8J^ 
h 



for < A < h< 2A 



(8.1) 
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and 



472 

\,h ~ for < 2A < h . (8.2) 



The value is just the activation energy I\^- for the transition between — 1 and 
0. The activation energy for the transition between and +1 is always (approximately) 
given by: 

In region III we have r A ^- > r^ + - and this is the reason for (8.2); but in Region IV 
we have the opposite I\^- < r^ + - so that we get : 

4J2 

Tx,h ~ j— j for < -A < h . (8.4) 

This answers a question raised in [R] about the "validity of Van't Hoff-Arrhenius law" 
which would predict, in our case a decay — 1 — > — > +1 with an asymptotics of the 
transition time determined by T~-~^-. 

Our results can be interpreted by saying that this law is valid in the region III whereas 
it is violated in regions II and IV for different reasons. 

The new phenomenon about which apparently there is no reference even in the physics 
literature is the possibility of a "direct" transition between —1 and +1 and also a possible 
change in the mechanism of transition for different values of the parameters. 

Notice that if we take a fixed small value of h and we vary A the analytic expression of 
changes when we cross the lines h = 2 A, A = 0. 

We draw in Fig. 8.1 the graph of r^,^ as a functin of A for a fixed small value of h. 
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Appendix. 

In this appendix we want to state and prove Proposition Al below. It refers to the 
first escape from a transient cycle A (see below) and, roughly speaking, it says that, under 
general hypotheses, with high probability, after many attempts, soon or later our process 
will really escape from A entering into another different cycle by passing through one of 
the minimal saddles of the boundary of A. 

The time for this transition has about the same asymptotics as the first hitting time to 
the boundary of A. 

The result of Proposition Al was already used before without an explicit proof (see, 
e.g., [KOI] and [K02]); it is, in fact, a simple consequence of the strong Markov property 
but we think it useful, in order to better explain its statement, to eventually provide an 
explicit proof. 

We will state our results in a slightly more general set-up than the one considered in 
the present work (we also use a different notation). We will consider general Metropolis 
reversible Markov chains. 

We suppose given an ergodic, aperiodic Markov chain (^)t=o,i,2,... with finite state- 
space O and with transition probabilities P(x,y) satisfying the following 

Hypothesis M. 

There exists a function H : fi — > R + such that 

P(x, y) = q(x, y) exp(-/3[H(y) - H(x)}+) , (Al) 
where q(x, y) = q(y, x) and (a)+ is the positive part (:= max{a, 0}) of the real number a. 



The above choice corresponds to a Metropolis Markov chain which is reversible in the 
sense that: 

\fx,x' GO : /i(x)P(x,x') = n(x')P(x',x) (A2) 

with 

H(x) oc exp(-(3H(x)) . (A3) 

One can introduce the notions of pair of communicating states, path, connected subset of 
O, boundary dQ of set Q C O cycles, ... as the obvious generalizations of the corresponding 
ones given in Section 2. 

For any set Q C O we introduce the set of all the minima of the energy in the boundary 
dQ of Q: 

U(Q) := {zedQ: min H{x) = H(z)} . (A4) 

xEdQ 
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By F(Q) we denote the set of the absolute minima of the energy in the set Q C 0: 

F(Q) :={yeQ: mmH(x) = H(y)} . (A5) 

xeQ 

A cycle A for which there exists y* G U(A) downhill connected to some point x in A c 
(namely 3 x G" A, communicating with y*, with H{x) < H(y*) =: H(U(A))) : is called 
transient; points like y* are called (minimal) saddles. S(A) will denote the set of all 
minimal saddles of A. 

Let R = R(A) be the subset of A to which some point in S(A) is downhill connected: 

R(A) := {y eA such that 3z G S(A) with P(x, y) > 0} ; (A6) 

let V = V(A) be the analogue of R outside A: 

V(A) :={y<£A such that 3z E S(A) with P(x, y) > 0} . (A 7) 

We set: 

H := R(A) U V(A) . (A.8) 

Proposition A.l. 

Consider a transient cycle A. Given e > let 

T(e):=expp[H(S(A))-H(F(A)) + e] . (A. 9) 

Then, for every e > 0, i6A, 

lim P x {T {A[JdA y > T{e)) = , (A10) 



and 



Proof. 



hm P,(X T(AuaA)c G = 1 . (All) 



From Hypothesis M and the definition of 5(A) we know that there exists a positive 
constant c > 0, independent of /?, such that 

inf P(x,y) > c, lim sup P(x 7 y) = . (A12) 

a;e5(A),ye-H ' P^°° x€S(A),y2H 
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We define, now, the sequence T\ of stopping times corresponding to subsequent passages 
of our chain X t in dA: 

t :=m£{t >0:X t edA} 
en :=m£{t > t :X t £ dA} ' 



and for j = 1,2,.. .: 



(A13) 



(A14) 



Tj :=inf{t > a i : X t G dA} 
cfj :=inf{t > Tj-i : X t g dA} 

We set, for j = 1,2,. . .: 

z 3 = fo-i + 1, t j] , T j : = \ = t j ~ ~ 1 • ( Al5 ) 
Suppose that X Tj _ 1+ i G A; let 

a* := min{t > r, : X t ^ X T .} ; (A16) 

we say that the interval X, is good if the following conditions are satisfied: 

T S < T(e) 
X T . G S(A) . 
X a * G H 

3 

Let 

j* := min{j : Tj is not good} . 

Given the integer N we want to estimate, for every x G A, the probability P x (tv(A) > 
NT(e)). 
We write: 

Px(t V (a) > NT(e)) = P x (rv(A) > NT(e) ; j* > N) + P x (t v{a) > NT(e) ; j* < N) 

(A17) 

Let us consider the first event in the decomposition given in (A. 17): {t V (a) > NT(e); j* > 
N}. 

We have: 

P x (t v(a) > NT(e) ; j* > N) < 
<P X (X T1 G S(A); X T1+1 ?V(A),...,X TN G S(A); X TN+1 ?V(A)) < 
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< (l-c) N . (A.18) 

Now, from Proposition 3.7 of [OS1], (A. 12) and the stationarity of our Markov chain we 
have : 

P x (Xj is not good) < 6(J3) (A. 19) 



with 



lim 6(/3) = . (A.20) 



From (A.20) we get: 

PArv(A) > NT(e) ; j* < N) < P x (j* < N) < 6(0) N . (A.21) 
To conclude the proof it suffices to choose : 

N = N((3) = l/6((3). 
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